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STUDIES  IN  GRAVIMETRIC  GEODESY 

1.  Introduction 

This  is  a final  report  and  summary  of  the  work  done  under 
the  Contract  No.  F 19628 - 77-C -0082.  The  contractual  and  reporting 
period  is  1 January,  1977  - 30  September,  1978.  Objectives  of 
the  work  as  defined  in  the  contract  were  as  follows: 

a.  Determine  where  a limited  number  of  new  absolute 
gravity  measurements  will  have  the  largest  influ- 
ence on  improving  the  IGSN  7 1 network;  determine 
proper  intervals  for  additional  absolute  sites  in  order 
to  solve  the  second  and  higher  order  correction  terms 
for  calibration  of  gravimeters;  determine  the  locations 
for  new  absolute  measurements  which  reduce  systematic 
errors  to  a minimum. 

b.  Study  the  Geodetic  Boundary  Value  problem  in  light 
of  the  theoretical  breakthrough  of  Lars  Hormander's 
"The  Boundary- value  Problem  of  Physical  Geodesy", 
published  in  Stockholm,  Sweden  in  1975. 

c.  Study  the  homogeneity  and  isotropy  assumptions  in 
gravity  predictions.  Form  models  that  vary  from  one 
area  to  another,  and  models  that  may  have  azimuth 
dependence  in  their  covariance  functions. 

The  work  done  on  items  b and  c have  been  reported  as 
Scientific  Reports  No.  2 and  No.  1 respectively.  The  work  done 
under  item  a has  been  reported  in  the  quarterly  status  reports  and 
in  informal  communications  sent  to  the  Monitor  of  the  Contract. 
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In  the  following  only  a short  summary  will  be  given  of  work 
done  under  items  b and  c because  the  detailed  technical  reports, 
submitted  ea  rlier,  cover  the  items  in  detail.  More  detailed  report 
will  be  presented  about  the  work  done  under  item  a. 

2.  Geodetic  bo unda  ry  value  problem 

The  work  done  in  this  problem  was  reported  in  the  techni- 
cal report  by  Helmut  Moritz  (1977):  "Recent  Developments  in  the 
Geodetic  Bounda  ry- Value  Problem."  This  report  was  submitted 
as  Scientific  Report  No.  2 under  this  contract  and  was  distributed 
by  the  Air  Force  Geophysic  Laboratory  as  document  AFGL-TR-78- 
0002.  It  was  also  distributed  as  Report  No.  266  in  the  series  of 
Reports  of  the  Department  of  Geodetic  Science. 

The  abstract  of  the  report  reads: 

"The  report  reviews  progress  in  the  mathematical  formula- 
tion and  treatment  of  the  geodetic  boundary- value  problem,  in 
particular,  the  existence  and  uniqueness  theorems  of  L.  Hormander 
and  the  gravity  space  approach  due  to  F.  Sanso.  The  method  of 
Hormander  uses  a very  advanced  inverse  function  theorem  of  non- 
linear functional  analysis.  Sanso  has  transformed  Molodensky's 
free  boundary- value  problem  into  a fixed  boundary- value  problem 
in  'gravity  space'  , thereby  essentially  reducing  the  mathematical 
complexity.  As  a linear  approximation,  the  gravity  space  approach 
gives  identical  results  to  the  conventional  linearization,  but  gravity 
space  appears  superior  for  treating  questions  of  existence  and 
uniqueness  of  the  solution,  although  it  is  restricted  to  the  pure 
gravitational  case  without  centrifugal  force." 

It  is  recommended  that  those  who  are  interested  in  this 
subject  matter  read  the  original  excellent  technical  report.  The 
importance  of  this  work  might  be  reflected  by  the  following  quota- 
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tion  from  Moritz's  report:  "The  impact  of  the  gravity  space  approach 
to  the  theory  of  Molodensky's  problem  appears  to  be  enormous;  it 
may  well  be  comparable  to  the  impact  of  Hamiltonian  methods  to 
Newtonian  classical  mechanics  (both  apply  a Legendre  transforma- 
tion) ).  " 

3.  Homogeneity  and  isotropy  assumptions  in  gravity  prediction 

The  work  done  in  this  problem  was  reported  in  the  technical 
report  by  William  Kearsley  (1977):  "Non-Stationa ry  Estimation  in 
Gravity  Prediction  Problems.  " This  report  was  submitted  as 
Scientific  Report  No.  1 under  this  contract  and  was  di  stributed  by 
the  Air  Force  Geophysics  Laboratory  as  document  AFGL-TR-77- 
0 186.  It  was  also  distributed  as  Report  No.  256  in  the  series 
Reports  of  the  Department  of  Geodetic  Science. 

The  abstract  of  the  report  reads: 

"This  report  investigates  the  impact  that  the  assumptions 
of  homogeneity  and  isotropy,  when  applied  to  potential  related 
fields,  have  upon  the  stochastic  processes  which  are  applied  to 
these  fields.  After  seeing  how  these  assumptions  are  incorporated 
into  the  statistical  model  to  produce  the  familiar  covariance  function, 
the  investigation  centers  on  techniques  which  can  be  used  to  detect 
the  presence  of  anisotropy  in  the  field.  The  method  found  most 
useful  in  the  two-dimensional  covariance  function,  and  some  methods 
of  representing  this  function  are  also  investigated. 

Numerical  studies  are  then  carried  out  to  see  the  effect  the 
use  of  the  2-D  covariance  function  has  upon  the  results  of  prediction 
and  collocation  computations.  It  is  found  that,  under  certain  circum- 
stances, the  2-D  function  produces  a result  superior  to  that  given 
by  the  general  function.  Recommendations  are  then  given  as  to 
when  the  2-D  covariance  function  should  be  used  in  practical  solu- 
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tions,  and  suggestions  made  as  to  the  possible  areas  of  further 
resea  rch. 1 ' 

More  details  are  reflected  in  the  section  of  Conclusions  in 
the  report  which  reads: 

"The  two-dimensional  covariance  function  provides  the  most 
efficient  means  of  detecting  and  representing  the  anisotropic  char- 
acteristics of  a data  set  distributed  over  a plane.  This  function 
graphically  describes  the  covariance  which  exi sts  between  point 
pairs  of  all  separations  and  orientations.  The  extent  of  anisotropy 
is  indicated  by  the  departure  of  the  contours  of  the  2-D  covariance 
surface  from  a circular  pattern,  and  the  orientation  of  the  axes  of 
maximum  and  minimum  correlation  are  clearly  shown. 

It  is  possible  to  model  the  2-D  covariance  surface  by  gen- 
erating a simple  covariance  function  for  each  azimuth,  0 s a < 360. 
The  logarithmic  function  suggested  by  Moritz  (1976,  p.  29)  appears 
to  be  the  best  model  overall,  particularly  when  the  function  attains 
negative  values.  The  fact  that  the  2-D  cross-covariance  function 
is  not  symmetrical  complicates  the  generation  of  the  surface  by 
this  method.  It  is  possible  to  overcome  this  problem  by  using  the 
symmetrical  2-D  function  to  approximate  the  cross-covariance  sur- 
face. 

The  ideal  function  would  enable  the  generation  of  all  auto 
and  cross-covariance  functions  knowing  the  pertinent  parameters  of 
(say)  the  anomalous  gravity  field.  The  third-order  Markov  function 
suggested  by  Jordan  (1972)  has  this  capability.  Unfortunately,  the 
theoretical  relationships  did  not  agree  with  the  actual  relationships 
in  this  instance.  It  is  felt  that  this  is  an  important  area  for 
further  research,  if  the  usefulness  of  the  2-D  covariance  function 
is  to  be  fully  exploited. 

The  2-D  covariance  function  is  capable  of  producing  results 
superior  to  those  obtained  by  the  general  function  when  ce:  tain 


4 


conditions  are  present.  These  conditions  will  produce  large  differ- 
ences between  elements  of  the  covariance  matrices  derived  from 
the  general  covardQce  analysis  and  from  the  2-D  covariance  analy- 
sis. They  will  occur  when: 

(i)  anisotropic  effects  are  present  and,  because  of  the  dis- 
tribution of  the  data,  predictions  must  be  performed 
over  large  separations  and  in  an  asymmetric  configur- 
ation, or 

(ii)  anisotropy  is  strongly  evident  and  homogeneous  through- 
out the  field.  Such  an  effect  can  be  seen  in  areas  where 
geoidal  slope  are  uniformly  and  consistently  large  (e.  g. 
the  geoid  slope  across  Australia).  In  fact,  under  these 
conditions  the  solution  using  the  general  function  appears 
to  break  down. 

In  any  case,  a 2-D  covariance  analysis  should  be  performed 
on  data  which  shows  anisotropic  tendencies.  This  will  indicate  the 
extent  to  the  azimuth  dependence  of  the  covariance  function  and 
enable  remedial  action  to  be  taken  (e.  g.  in  the  configuration  of  the 
data  used  in  subsequent  computation)  if  this  appears  warranted. 

The  2-D  covariance  surface  may  also  provide  useful  infor- 
mation concerning  a suitable  "trend  surface"  to  be  fitted  to  the 
original  data.  Knowing  that  the  residuals  of  the  actual  data  from 
the  trend  should  be  isotropic,  it  should  be  possible  to  discover 
what  nature  of  surface  must  be  fitted  in  order  to  transform  the  2-D 
covariance  surface  to  a surface  of  revolution.  (This  may  be  best 
performed  in  the  spectral  domain.)  The  residuals  can  then  be  used 
in  the  stochastic  processes  with  the  knowledge  that  they  do  in 
reality  possess  isotropic  characteristics." 

For  more  details  the  reader  is  referred  to  the  original 

report. 
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4.  Selection  of  preferred  locations  of  new  absolute  gravity 
measurements  in  the  gravity  networks 

4.  1 Introduction 

International  Gravity  Standardization  Net  197  1 [iGSN  7 l]  was 
adopted  by  the  International  Union  of  Geodesy  and  Geophysics  at 
the  XV  General  Assembly  in  Moscow  197  1 (Morelli,  et  al  1974). 

In  the  least  squares  adjustment  which  produced  the  IGSN  7 1 ten 
absolute  measurements  of  gravity  were  used  at  the  eight  sites. 

Since  the  adoption  of  IGSN  7 1,  new,  very  accurate,  portable  absolute 
gravity  measuring  devices  have  been  developed.  During  the  develop- 
ment phase  of  these  apparatuses,  a question  "where  is  the  best 
place  to  make  an  additional  absolute  gravity  determination  to  improve 
most  of  the  IGSN  71"  was  posed.  Under  this  research  project  the 
analyses  were  performed  in  order  to  answer  this  question. 

4.2  Mathematical  model 

The  IGSN  7 1 is  formed  by  1854  gravity  stations,  distributed 
around  the  world.  The  inverse  of  the  normal  matrix  for  the 
solution  of  the  gravity  values  of  the  stations  and  other  parameters 
was  available  to  us  and  obtained  from  the  Defense  Mapping  Agency, 

Aerospace  Center,  Geodetic  Survey  Squadron,  F.E.  Warren  AFB, 

Wyoming,  where  the  final  simultaneous  adjustment  of  IGSN  71  was 
made.  The  analyses  of  the  effect  of  new  absolute  measurements  to 
the  variance -co variance  matrix  can  be  accomplished  using  step-by- 
step  sequential  solutions.  A brief  outline  of  a procedure  is  given 
here  following  the  notations  and  derivations  given  by  Uotila  (1973a). 

Assuming  we  have  a set  of  observations  L*  which  are  functions 
of  a set  of  parameters,  Xa: 


L‘a  = F1(XJ. 


(1) 


xrx« 

where  X0  = approximate  values  of  the  parameters,  we  get  linearized 
observation  equations 

V,  = A,Xj  + Llt  (2) 


where  Vj  - residuals 


X,  - X.1  - X, 


X1  - adjusted  values  of  parameters 


Li  Ln  - L* 


- observed  values 


L‘  = F i(X0). 

The  minimum  variance  solution  gives  us: 

Xj  = -(A1!  P1Ai)'1  AT,  Pt  Lt  (3) 

where  P{'  = y j variance-covariance  matrix  of  observations 

Lb 


The  adjusted  values  of  parameters  are: 

X‘  = x0  + xt 


and  their  variance-covariance  matrix 


“ (-^Ti  Pi  Aj)  1 - Ni1 


If  we  have  the  second  set  of  observations  and  their  var- 
iance-covariance matrix  £ 2 = P£l  and  a mathematical  model: 

Lb 

L f = F(X  ) (5) 


■ 


then  a combined  solution  of  the  pa  rameters  is: 


Xf  = X0  + X2 


where 


X,  = -(A\  P,  A,  + AS  P2  A2)-‘ 
(AtjP,  A,  + At2  P2  L2) 


and 


A 

a SX 


X =Xn 


» L2  - Lq  - L.(j 


Lo  = F2(X0) 

The  variance-covariance  matrix  of  X2  is: 


Sxz  = (N,  + A^2P2  A2)-1 


(6) 


(?) 


(8) 


Equation  (8) 


or 


can  be  easily  modified  to  (Uotila,  1973b): 

= Nf‘  - Nj1  At2(A2  At2  + P2-‘)-‘  A2N,'‘  (9) 

a 

Fxi  - At2(A2  N il  A\  + P2l)'‘  A2N f1  (10) 


a a 


Equation  (10)  gives  the  difference  of  the  variance-covariance 
matrix  of  parameters  as  obtained  from  the  first  set  of  observations 
and  the  variance-covariance  matrix  of  parameters  obtained  using  the 
first  set  of  observations  and  the  second  set  of  observations  combined. 

In  our  analyses  could  be  the  variance-covariance  matrix  of  gravity 

a 

values  of  IGSN  7 1 stations  and  £^2  the  new  variance-covariance  matrix 
of  the  same  gravity  values  after  new  absolute  measurements  have  been 
added  to  the  net.  We  are  interested  in  this  change  in  order  to  make  the 
best  site  selections  for  new  absolute  measurements.  Obviously  the  same 
kind  of  sequential  solution  can  be  continued  by  adding  and  getting 
£X3  using  £x2  as  starting  matrix  and  so  on.  We  co  dd  also  add  one 
observation  at  a time.  For  example,  we  can  add  a new  observation  at 


I 
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each  station,  one  at  a time,  and  see  the  influence  of  each  added 
observation  to  the  original  y ..  If  we  take  all  1854  stations  we  would 

a 

have  equal  number  of  y . - £ 2 differences.  Now  we  must  make  a 

A A 

. a a 

decision,  which  one  of  the  new  observations  is  giving  the  optimum 

change  in  the  original  variance-covariance  matrix.  When  the 

selection  is  made  we  have  a new  y 2 for  the  net  including  a new 

X. 

a 

selected  absolute  measurement.  We  can  then  find,  using  similar 
techniques,  where  the  second  absolute  measurement  should  be  made 
using  the  difference  £^  - £^3.  But  how  do  we  select  the  optimum 

a a 

£ i - £ 2 from  the  1854  possibilities? 

a ^a 

According  to  Fedorov  (1972)  there  are  several  properties 
which  could  be  used  to  determine  which  one  of  the  two  experiments 
is  the  preferred  one: 

a)  Experiment  E[  is  preferred  to  experiment  E2  if  the 

difference  £ - £ is  a positive-definite  matrix  or  in  other  words, 

Ei  > E,  if  £ < £ where  £ and  £ are  variance-covariance 

Ei  Js2  E3  E2 

matrices  of  the  corresponding  results  of  the  experiments. 

b)  The  second  criterion  is 

E‘  > E*  if  I^Ej  I < 

where  jv  I and  |£  I are  determints  of  the  variance-covariance 
Ei  E2 


matric  e s. 

c)  The  third  criterion  is: 


Et  > E2  if  Tr  £ < Tr  £ 

Ei  E2 


where  Tr  £ and  Tr  £ are  traces  of  the  variance-covariance 
Ei  E2 

matrices. 


d)  The  fourth  criterion  is: 

Ei  > E2  if  the  maximum  variance  of  E3  < maximum  variance  of  E2. 
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e)  E(  > E2  if  the  variance  of  a function  of  Et  < the  variance 
of  the  same  function  of  E2.  We  could  continue  with  this  criteria 
of  the  function,  including  a - d criteria  to  a set  of  functions,  but 
the  above  is  sufficient  to  show  that  there  are  several  possibilities 
for  selecting  criteria  to  decide  which  experiment  is  the  "optimum." 

Some  test  analyses  were  done  and  it  was  found  that  b and  c 
of  the  criteria  were  not  giving  much  different  selections.  The 
criteria  c and  d are  fast  computationally,  but  d reflects  a local 
improvement  and  not  necessarily  global,  therefore  the  c-criterion 
was  selected  to  be  used  in  these  analyses. 

If  we  have  three  matrices  A,  B and  C of  the  same  order 

and 

A - B = C 

then 

Tr(A)  - Tr(B)  = Tr(C). 

Letting  A = e j,  B = and  C = N,-‘At2  [A.N^  + P£1]'1  A2N{1,  it 
can  be  seen  that  the  smallest  Tr  (B)  is  produced  when  the  Tr  (C ) 
is  largest  since  Tr  (A)  is  invariant  in  this  case.  Thus,  the  pro- 
blem of  finding  the  optimum  site  for  the  first  new  absolute  measure- 
ment is  to  find  which  one  of  the  added  observations  maximizes  the 
trace  of  C. 

4.  3 Computational  technique 

If  a single,  uncorrelated  observation,  which  is  a direct  ob- 
servation of  a parameter,  is  added  to  the  system,  it  turns  out  that 
the  computation  of  the  change  in  the  trace  is  relatively  fa  st  and  easy. 
Let's  assume  that  a new  absolute  measurement  of  gravity  is  done 

at  station  i.  We  wish  to  evaluate  the  right  side  of  the  equation  (10). 

, , .th 

A2-matnx  is  a row  matrix  having  zero  elements  except  at  l 

column  there  is  +1,  therefore  A2N“lA^  is  a number  equal  to  the 
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station  and  P2* 


variance  of  gravity  value  of  i 


is  a number  equal 


to  the  variance  of  the  new  absolute  measurement  at  the  i station. 

Therefore,  A2N11AT2+  P2*  is  a sum  of  these  two  variances  and  its 

inverse  is  the  reciprocal  of  this  sum.  The  matrix  product  N,*AT 

in  this  case  will  be  equal  to  the  i*^  column  of  the  N, ''-matrix  and 
th 

A2N('  will  be  the  i row  of  the  Nj1  matrix;  therefore,  for  tlu 


change  of  the  trace  of  C matrix  (th<  right  side  of  equation  10 i,  we 

need  the  sum  of  the  squares  of  the  elements  of  i^1  row  of  |! 

matrix  multiplied  by  the  reciprocal  of  sum  of  the  variant  < of  the  ra\  tv 
i c th 

value  of  1 station  and  the  variance  of  the  new  absolute  measurement. 

The  computational  procedure  described  above  turned  >>ut  t.  he 
simple;  therefore  it  is  feasible  to  compute  the  change  in  the  trac  <• 
corresponding  for  each  case,  where  each  one  of  the  stations  was 
occupied  by  an  ab  solute  appa  ratus  and  a new  absolute  measurement 
of  gravity  was  performed.  The  "optimum"  station  to  be  occupied 
first  would  be  the  station,  which  produced  a maximum  change  in 
the  trace  as  described  above. 


4.4  Variance-covariance  matrix  of  1GSN  7 1 gravity  values 


The  inverse  of  the  normal  matrix  received  from  Geodetic  Survey 
Squadron  was  not  the  variance-covariance  matrix  of  gravity  values 
of  IGSN  7 1 stations,  but  a weight  coefficient  matrix.  There  was  an 
unknown  scale  factor  involved.  It  was  solved  by  comparing  variances 
derived  from  standard  deviations  given  in  (Morelli,  et  al,  1974)  and 
diagonal  elements  of  the  matrix.  Comparing  283  selec  ted  va  rianc e s 
in  IGSN  7 1 with  corresponding  elements  of  the  matrix,  it  was  found 
that  the  scale  factor  was  0.0037788.  After  multiplying  the  weight- 
coefficient  matrix  by  0.0037788,  we  obtained  a variance-covariance 
matrix  corresponding  to  the  standard  errors  given  in  (Morelli,  et 
al,  1974).  A brief  examination  of  standard  errors  given  for  IGSN 
7 1 stations  revealed  that  the  standard  errors  seem  to  be  too  small. 


For  example,  Washington,  11687  R ha  s a standard  deviation  0.0  11 
mgal.  The  major  contributor  to  the  standard  errors  is  the  limited 
number  of  absolute  measurements  and  thei  r di  stribution  m the  net. 

Even  if  all  ten  ab  solute  mea  surements  used  in  the  adjustment  of 
1GSN  7 1 were  done  at  a single  site  with  accuracies  stated  in  the 
publication  (Morelli,  et  al,  1974),  the  standard  error  of  the 
weighted  mean  would  have  been  0.0  143  mgal;  therefore  it  is 
impossible  to  have  a standard  error  less  than  this  value  for  any 
gravity  value  in  IGSN  7 1 net.  A closer  look  of  weighting  systems 
used  in  the  final  IGSN  7 1 suggested  that  there  might  have  been  a 
problem  in  relative  weighting  of  the  absolute  measurements  with  other 
measurements.  If  that  was  the  case,  the  gravity  values  of  IGSN  7 1 
would  not  change  but  variances  for  them  would  increase. 

We  had  available  a variance-covariance  matrix  of  gravity 
values  for  372  gravity  stations  computed  by  Uotila  using  linear 
correction  term  to  calibration  of  gravimeters  under  AFC  RL  Contract 
No.  F 19628-68-C-0335.  Comparing  variances  of  366  common  stations 
it  was  found  that  the  above  obtained  variance-covariance  matrix  of 
IGSN  7 1 should  be  multiplied  by  1.90848  in  order  to  have  a reasonable 
agreement  between  these  two  variance-covariance  matrices.  There- 
fore, the  original  weight  coefficient  matrix  of  IGSN  7 1 wa  s multiplied 
by  0.0072118  16  in  order  to  obtain  the  variance-covariance  matrix  of 
gravity  values  of  IGSN  7 1.  That  means  that  standard  deviations 
given  in  IGSN  7 1 publications  should  be  multiplied  by  >/l.  90848~  1.  38 , 
or,  in  other  words,  should  be  increased  38/°. 

4.  5 Selection  of  the  network  to  be  analyzed 

As  was  mentioned  above  there  are  1864  stations  in  IGSN  7 1 
network.  Many  of  these  stations  having  the  same  IC  B-identification 
number  are  so  called  excentric  stations  and  are  highly  correlated 
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with  one  another.  An  improvement  in  one  would  result  an  improve- 
ment in  all  excentric  stations  and  the  improvement  in  trace  would 
reflect  a strong  local  influence  rather  than  a global  one.  To 
reduce  this  local  dependency  on  the  trace,  a net  of  422  stations 
was  selected  from  the  1854  IGSN  7 1 net.  The  criterion  for  their 
selection  was  that  no  more  than  one  station  from  given  IG  B number 
was  included.  Within  a given  IGB  number  the  station  with  the 
most  external  ties  was  selected.  See  Table  1 listing  of  the  stations 
and  Figure  1 for  the  distribution  of  stations. 

4.6  Results  of  analyses 
4.6.  1 World  wide  analyses 

A few  test  runs  were  made  with  the  net  containing  these  422 
stations  and  assuming  n = 0.  02  mgal  accuracy  for  new  absolute 
measurements  of  gravity.  It  was  noticed  that  the  stations  having 
very  large  variance  were  selected  as  the  best  candidates  for  new 
absolute  measurements.  If  these  stations  were  really  selected, 
there  would  be  hardly  any  improvements  elsewhere  in  the  net 
except  in  the  variance  at  these  stations;  therefore,  the  procedure 
was  modified.  We  subtracted  from  the  trace  the  improvement  of 
the  station,  where  the  absolute  measurement  was  made.  This 
"partial  trace"  reflected  better,  in  our  opinion,  an  improvement 
in  the  whole  global  net. 

The  results  of  the  new  analyses  gave  the  following  priorities 
for  new  absolute  measurements  (using  0.02  mgal  accuracy  for  the 
new  absolute  measurements): 


1. 

Nairobi 

357  16 

A 

2. 

N ew  D el  hi 

10  187 

K 

3. 

Rio  Gallegos 

5 1119 

K 
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TABLE  1 


GRAVITY  BASE  STATIONS  USED  IN  THE  ANALYSIS 


U.S.  NETS 


List  Ing 

In 

IGB  NUMBER  order 

CODE 

IGB 

B141 

0 

KEY  WEST 

8150 

R 

MIAMI 

8160 

J 

WEST  PALM  BEACH 

8170 

K 

VERO  BEACH 

U 

8172 

J 

TAMPA 

8180 

J 

COCOA 

3181 

K 

ORLANDO 

I 

3191 

J 

DAYTONA  BEACH 

u 

819  1 

0 

DAYTONA  BEACH 

u 

8277 

J 

CORPUS  CHRIST 

8279 

J 

LAREDO 

I 

8289 

B 

COTULLA 

8290 

J 

NEW  ORLEANS 

8295 

J 

HOUSTON 

8298 

M 

SAN  ANTONIO 

u 

1 1629 

J 

CHARLESTON  J 

I 

1 1629 

L 

CHARLESTON 

1 1649 

J 

FI.ORENCE/S.  CAROLINA® 

1 1658 

J 

RALEIGH 

11677 

J 

RICHMOND 

1 1687 

M 

WASHINGTON 

1 1701 

J 

JACKSONVILLE 

I 

11711 

J 

BRUNSWICK 

u 

11711 

K 

BRUNSWICK 

1 1714 

J 

ALBANY 

1 1720 

.1 

BEAUFORT 

1 1721 

J 

SAVANNAH 

I 

1 1734 

J 

ATLANTA 

u 

1 1734 

K 

ATLANTA 

1 1750 

J 

CHARLOTTE 

1 1753 

J 

KNOXVILLE 

1 1759 

J 

MEMPHIS 

I 

11785 

J 

LOUISVILLE 

1 1807 

K 

AUSTIN 

1 1826 

J 

DALLAS 

1 1842 

J 

LITTLE  ROCK 

1 1877 

J 

WICHITA 

I 

1 1880 

L 

ST.  LOUIS 

u 

1 1880 

M 

ST.  LOUIS  M M 

I 

1 1894 

J 

KANSAS  CITY 

V 

1 1894 

K 

KANSAS  CITY 

1 19  16 

J 

EL  PASO 

I 

1 1926 

J 

ALAMOGORDO 

I 

1 1931 

J 

LUBBOCK 

1 195  1 

J 

AMARILLO 

1 1956 

J 

A I.  B UftUE  ROUE 

I 

1 1994 

A 

DENVER 

U 

1 1994 

N 

DENVER  N CO 

1 1998 

J 

GRAND  JUNCTION 

12027 

K 

SAN  DIECO 

12032 

J 

PHOENIX 

12038 

K 

LOS  ANGELES 

u 

12047 

K 

NORTON  AFB  K 

12065 

J 

LAS  VEGAS 

12099 

J 

RENO 

12172 

0 

SAN  FRANCISCO 

12181 

J 

FAIRFIELD 

15148 

J 

BANGOR 

15  167 

J 

CARIBOU 

15203 

R 

NEW  YORK  CITY 

15204 

J 

PRINCETON 

15209 

J 

PITTSBURG 

Listing  In  ALPHABETICAL  order 
IGB 

CODE  NUMBER  NAME 


I 

I 1926 

J 

ALAMOGORDO 

1 1714 

J 

ALBANY 

1 1956 

J 

ALBUQUERQUE 

11951 

J 

AMARILLO 

I 

1 1734 

J 

ATLANTA 

U 

1 1734 

K 

ATLANTA 

1 1807 

K 

AUSTIN 

15148 

J 

BANGOR 

1 1720 

J 

BEAUFORT 

I 

15558 

L 

BILLINGS 

u 

15558 

M 

BILLINGS  M MT 

15560 

K 

BISMARCK 

15636 

J 

BOISE 

15221 

J 

BOSTON 

I 

11711 

J 

BRUNSWICK 

u 

11711 

K 

BRUNSWICK 

15228 

J 

BUFFALO 

15167 

J 

CARIBOU 

15526 

L 

CASPER 

I 

1 1629 

L 

CHARLESTON 

u 

1 1629 

J 

CHARLESTON  J 

11750 

J 

CHARLOTTE 

15514 

M 

CHEYENNE 

15317 

M 

CHICAGO 

8180 

J 

COCOA 

15303 

J 

COLUMBUS 

u 

8277 

J 

CORPUS  CHRIST 

I 

8289 

B 

COTULLA 

u 

15682 

B 

CUTBANK  B 

1 1826 

J 

DALLAS 

I 

8191 

J 

DAYTONA  BEACH 

u 

8191 

0 

DAYTONA  BEACH 

I 

1 1994 

A 

DENVER 

u 

1 1994 

N 

DENVER  N CO 

15323 

J 

DETROIT 

I 

15462 

J 

DULUTH 

1 19  16 

J 

EL  PASO 

12181 

J 

FAIRFIELD 

I 

15466 

J 

FARGO 

1 1649 

J 

FLORENCE/'S . CAROL  1 NA® 

15416 

J 

FREMONT 

I 

15477 

M 

GRAND  FORKS 

1 1998 

J 

GRAND  JUNCTION 

15671 

L 

GREAT  FALLS 

8295 

J 

HOUSTON 

1 1701 

J 

JACKSONVILLE 

I 

1 1894 

J 

KANSAS  CITY 

u 

1 1894 

K 

KANSAS  CITY 

8141 

O 

KEY  WEST 

1 1753 

J 

KNOXVILLE 

8279 

J 

LAREDO 

12065 

J 

LAS  VEGAS 

1 1842 

J 

LITTLE  ROCK 

12038 

K 

LOS  ANGELES 

I 

1 1785 

J 

LOUISVILLE 

I 

1 1931 

J 

LUBBOCK 

I 

15339 

A 

MADISON 

u 

15339 

J 

MADISON  J WI 

I 

15722 

J 

MEDFORD 

1 1759 

J 

MEMPHIS 

8150 

R 

MIAMI 

I 

15212 

J 

MIDDLETOWN 

CODE 


I = > IGSN71  net  only,  V ->  VAV  net  only,  NO  CODE  =>  common  to  both  nets 


TABLE  ! (continued) 

U.S.  NETS  (continued) 


Lis  t ing 

i n 

IGB  NUMBER  order 

Lis  t Ing 

1 n 

ALPHABETICAL  order 

1GB 

IGB 

U 

15212 

A 

MIDDLETOWN  A 

U 

15212 

A 

MIDDLETOWN  A 

I 

15212 

J 

MIDDLETOWN 

15443 

L 

MINNEAPOLIS 

15221 

J 

BOSTON 

I 

15581 

L 

MINOT 

15220 

J 

BUFFALO 

8290 

J 

NEW  ORLEANS 

15230 

J 

PORTLAND/ MAINE® 

15203 

R 

NEW  YORK  CITY 

15236 

J 

SYRACUSE 

U 

12047 

K 

NORTON  AFB  K 

15303 

J 

COLUMBUS 

I 

156  1 1 

J 

OGDEN 

15317 

M 

CHICAGO 

8181 

K 

ORLANDO 

15323 

J 

DETROIT 

12032 

J 

PHOENIX 

I 

15339 

A 

MADISON 

15209 

J 

PITTSBURG 

u 

15339 

J 

MADISON  J XI 

15230 

J 

PORTLAND/  MAINE® 

u 

15414 

J 

STUART 

15752 

J 

PORTLAN  D/OREGON® 

15416 

J 

FREMONT 

15204 

J 

PRINCETON 

15426 

J 

SIOUX  CITY 

1 1658 

J 

RALEIGH 

15436 

.1 

SIOUX  FALLS 

15543 

J 

RAPID  CITY 

15443 

L 

MINNEAPOLIS 

12099 

J 

RENO 

I 

15462 

J 

DULUTH 

1 1677 

J 

RICHMOND 

I 

15466 

J 

FARGO 

U 

15601 

K 

SALT  LAKE  CIT 

I 

15477 

M 

GRAND  FORKS 

I 

1 560  I 

J 

SALT  LAKE  CITY 

15514 

M 

CHEYENNE 

8298 

N 

SAN  ANTONIO 

15526 

L 

CASPF.R 

12027 

K 

SAN  DIEGO 

15543 

J 

RAPID  CITY 

12172 

0 

SAN  FRANCISCO 

15546 

J 

SHERIDAN 

1 172  1 

J 

SAVANNAH 

I 

15550 

L 

BILLINGS 

15772 

P 

SEATTLE 

u 

15550 

M 

BILLINGS  M MT 

15546 

J 

SHERIDAN 

15560 

K 

BISMARCK 

15426 

J 

SIOUX  CITY 

1 

15581 

L 

MINOT 

15436 

J 

SIOUX  FALLS 

I 

15601 

J 

SALT  LAKE  CITY 

I 

15677 

J 

SPOKANF. 

u 

15601 

K 

SALT  LAKE  CIT 

I 

1 1 880 

L 

ST.  LOUIS 

I 

1561  1 

J 

OGDEN 

U 

1 1 880 

M 

ST.  LOUIS  M M 

15636 

J 

BOISE 

u 

15414 

J 

STUART 

15671 

L 

GREAT  FALLS 

15236 

J 

SYRACUSE 

I 

15677 

J 

SPOKANE 

u 

8172 

J 

TAMPA 

u 

15682 

B 

CUTBANK  B 

8170 

K 

VERO  BEACH 

I 

15722 

J 

MEDFORD 

1 1687 

M 

WASHINGTON 

15752 

J 

PORTLAND/OREGON® 

8160 

J 

WEST  PALM  BEACH 

15772 

P 

SEATTLE 

1 1877 

J 

WICHITA 

WORTH  AMERICAN  NETS  excluding  U.S, 

NETS 

U 

889 

A 

PANAMA  A 

4669 

K 

ACAPULCO 

I 

089 

M 

PANAMA 

I 

231  19 

K 

ANCHORAGE 

U 

899 

J 

CRISTOBAL 

U 

231  19 

J 

ANCHORAGE  J 

994 

K 

SAN  JOSE 

26703 

J 

BARTER  ISLAND 

4526 

K 

MANAGUA 

19214 

J 

CALGARY 

4539 

K 

SAN  SALVADOR 

2236  1 

J 

CAPE  DYER 

4640 

K 

GUATEMALA 

u 

19084 

R 

CHURCHILL  R 

U 

4640 

n 

GUATEMALA  M 

u 

899 

J 

CRISTOBAL 

4669 

K 

ACAPULCO 

u 

23049 

K 

DAWSON 

4698 

A 

PASO  DE  CORTES 

19233 

M 

EDMONTON 

4699 

A 

MEXICO  CITY 

26195 

K 

EUREKA 

8320 

K 

SAN  LUIS  POTOSI 

23147 

K 

FAIRBANKS 

8350 

K 

MONTERREY 

18788 

J 

FORT  CHI  MO 

15239 

J 

TORONTO 

u 

19382 

J 

FORT  NELSON 

15253 

J 

MONTREAL 

19360 

L 

FORT  ST.  JOHN 

15255 

J 

OTTAWA 

22338 

J 

FROBISHER  BAY 

15261 

J 

QUEBEC 

18730 

J 

GOOSE  BAY 

15282 

J 

ROBF.RVAL 

19258 

K 

CRANDF.  PRAIRIE 

15497 

0 

WINNIPEG 

4640 

K 

GUATEMALA 

15692 

A 

LETHBRIDGE 

u 

4640 

M 

GUATEMALA  M 

0 

15793 

J 

VANCOUVER  J 

22581 

J 

HALL  BEACH 

f 

15793 

M 

VANCOUVER 

15692 

A 

LF.THBR  1 DGF. 

18730 

J 

GOOSE  BAY 

I 

22485 

J 

LONGSTAFF 

CODE: 

I => 

' ICSN71  net  only.  U => 

UAU  net  onl 

ly.  NO 

CODE  = > common  to  bo 
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TABLE  1 (continued) 


NORTH  AMERICAN  NETS  excluding  U.S.  NETS  (continued) 

Listing  in  ICB  NUMBER  order  Listing  in  ALPHABETICAL  order 


IGB 

CODE  NUMBER  NAME 


IGB 

CODE  NUMBER  NAME 


1B746 

J 

SCHEFFERV I LLE 

4526 

K 

MANAGUA 

187  Bit 

J 

FORT  CHI MO 

4699 

A 

MEXICO  CITY 

u 

19084 

R 

CHURCHILL  H 

8350 

K 

MONTERREY 

19214 

J 

CAI.GARV 

15253 

J 

MONTREAL 

19223 

A 

RED  DEER 

26469 

0 

MOULD  BAY 

19233 

M 

EDMONTON 

U 

23121 

B 

NORTHWAY 

19258 

K 

GRANDE  PRAIRIE 

15255 

J 

OTTAWA 

19360 

1. 

FOR’r  ST.  JOHN 

I 

889 

M 

PANAMA 

V 

1 9382 

J 

FORT  NELSON 

U 

889 

A 

PANAMA  A 

22338 

J 

FROBISHER  BAY 

4698 

A 

PASO  DE  CORTES 

2236  1 

J 

CAPE  DYER 

26816 

A 

POINT  BARROW 

I 

22485 

J 

LONGSTAFF 

15261 

J 

QUEBEC 

22581 

J 

HALL  BEACH 

19223 

A 

RED  DEER 

22908 

L 

WATSON  LAKE 

26244 

K 

RESOLUTE  BAY 

23005 

A 

WHITEHORSE 

15282 

J 

ROBERVAL 

V 

23049 

K 

DAWSON 

994 

K 

SAN  JOSE 

V 

231  19 

J 

ANCHORACE  J 

8320 

K 

SAN  LUIS  POTOSI 

I 

231  19 

K 

ANCHORAGE 

4539 

K 

SAN  SALVADOR 

u 

23  1 20 

A 

SNAG 

18746 

J 

SCHEFFERV I LLE 

u 

23121 

B 

NORTHWAY 

U 

23120 

A 

SNAG 

23147 

K 

FA  I RBANKS 

15239 

J 

TORONTO 

26  195 

K 

EUREKA 

I 

15793 

M 

VANCOUVER 

26244 

K 

RESOLUTE  BAY 

u 

15793 

J 

VANCOUVER  J 

26469 

0 

MOULD  BAY 

22908 

L 

WATSON  LAKE 

26703 

J 

BARTER  ISLAND 

23005 

A 

WHITEHORSE 

26816 

A 

POINT  BARROW 

15497 

0 

WINNIPEG 

WORLD  NETS 

excluding  NORTH  AMERICAN  NETS 

150 

K 

ACCRA 

u 

41752 

J 

MARYBOROUGH 

I 

154 

L 

ABIDJAN 

I 

154 

L 

ABIDJAN 

I 

260 

K 

MONROVIA 

150 

K 

ACCRA 

I 

293 

.1 

CONAKRY 

3398 

K 

ADDIS  ABABA 

655 

J 

PARAMARIBO 

6824 

J 

ADEN 

668 

J 

GEORGETOWN 

10909 

J 

ACAD  I R 

U 

793 

K 

MATURIN 

18040 

J 

ACEN 

826 

K 

POPAYAN 

10177 

J 

AGRA 

836 

K 

CALI 

10132 

J 

AH  ME  DAB  AD 

844 

K 

BOGOTA 

I 

45466 

K 

ALB1JRY 

865 

K 

MEDELLIN 

I 

29522 

J 

ALERT 

U 

2087 

J 

KWAJALEIN 

I 

14463 

J 

ALGIERS 

2613 

A 

SINGAPORE 

I 

' *385 

T 

ALI  TERME 

U 

2622 

J 

MALACCA 

7 

- .933 

J 

ALICE  SPRINGS 

263  1 

J 

KUALA  LUMPUR 

093 

J 

ALTA 

2650 

J 

PENANG 

I 

3714 

J 

AMRITSAR 

2670 

J 

SONCKHLA 

u 

, 7904 

J 

ANGRI  J 

2969 

J 

COLOMBO 

i 

14192 

J 

ANKARA 

3302 

J 

ENTEBBE 

4371 

J 

ANTIGUA 

3398 

K 

ADDIS  ABABA 

40430 

K 

ANTOFAGASTA 

I 

3548 

J 

BANGUI 

i 

21572 

K 

APELVIKSAAS 

3609 

B 

LIBREVILLE 

u 

21572 

J 

APELVIKSAAS  J 

3649 

J 

DOUAI.A 

3686  1 

K 

ARF.QUIPA 

3663 

J 

LAGOS’ 

36880 

K 

ARICA 

I 

3728 

J 

BAMAKO 

32674 

J 

ASCENSION  ISLAND 

3836 

J 

BATHURST 

I 

6958 

K 

ASMARA 

3846 

J 

MBOUR- DAKAR 

V 

6958 

A 

ASMARA  A 

3885 

J 

NOUAKCHOTT 

40257 

J 

ASUNCION 

U 

3962 

J 

GAPE  VERDE  IS 

I 

10542 

K 

ASWAN 

4301 

J 

PORT  OF  SPAIN 

u 

45164 

C 

AUCKLAND 

4306 

K 

CARACAS 

1 1 187 

J 

AZORES 

I 

4341 

J 

ST.  LUCIA 

I 

21510 

A 

BAD  HARZBURC 

4371 

J 

ANTIGUA 

u 

21510 

C 

BAD  HARZBURC 

I 

4374 

J 

ST.  CROIX 

I 

21609 

J 

BAD  HF.RSFELD 

CODE: 

1 =) 

ICSN71  net  only,  U = > UAU  net 

on 

ly,  NO 

CODE  =>  common  to  bo 
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TABLE  I (continued) 


WORLD  NETS  excluding  NORTH  AMERICAN  NETS  (continued) 


Listing  in  IGB  NUMBER  order  Listing  in  ALPHABETICAL  order 


CODE 

1GB 

IGB 

4386 

J 

SAN  JUAN 

21500 

L 

BAD  NEUSTADT 

4387 

K 

RAMEY 

43982 

K 

BAHIA  BLANCA 

4404 

.1 

BAHRANQUILLA 

I 

3728 

J 

BAMAKO 

1 

•476 

J 

K 1 NCSTON 

17990 

B 

BAMBERG  SUD 

U 

476 

1. 

KINGSTON  L 

10052 

J 

BANARAS 

I 

4482 

J 

fOlir  AU  PRINCE 

I 

6537 

J 

BANGALORE 

1 

4407 

J 

MONTEGO  DAY 

6230 

J 

BANGKOK 

1 

4495 

J 

GUANTANAMO 

I 

3548 

J 

BANGUI 

u 

5295 

J 

HAWAII  ISLAND 

U 

6230 

M 

BANKOK 

u 

5696 

J 

WAKE.  ISLAND  J 

18012 

J 

BARCELONA 

I 

5696 

N 

WAKE  ISLAND 

25198 

K 

BARDUEOSS 

5854 

N 

GUAM 

4404 

J 

BARRANQUILLA 

6050 

1. 

MANILA 

M836 

J 

BATHURST 

6206 

J 

SAIGON 

46622 

J 

BEAUFORT  WEST 

6230 

.1 

BANGKOK 

14135 

J 

BE  I RUT 

u 

6230 

M 

BANKOK 

329  18 

L 

BELEM 

6366 

K 

RANGOON 

I 

36593 

J 

BELO  HORIZONTE 

I 

6430 

J 

MADRAS 

1 1524 

J 

BERMUDA 

I 

6537 

J 

BANGALORE 

u 

1 1524 

K 

BERMUDA 

I 

6578 

J 

HYDERABAD 

I 

17905 

L 

BIVIO  GIUNGANO 

I 

6592 

J 

BOMBAY 

25  1 74 

J 

BODO 

6824 

J 

ADEN 

844 

K 

BOGOTA 

6952 

K 

KHARTOUM 

I 

6592 

J 

BOMBAY 

I 

6956 

J 

TESSENE I 

u 

21523 

W 

BRANDENBURG 

u 

6958 

A 

ASMARA  A 

21520 

C 

BRAUNSCHWE I G 

I 

6950 

K 

ASMARA 

I 

36557 

J 

BRAZ ILIA 

6997 

K 

PORT  SUDAN 

I 

21638 

K 

BREMEN 

I 

7228 

J 

KANO 

41773 

J 

BRISBANE 

I 

7232 

J 

N 1 AMF.Y 

46603 

J 

BRISTOVN 

7407 

J 

PORT  ETIENNE 

I 

21604 

S 

BRUSSELS 

7485 

J 

GRAND  CANARY 

u 

21604 

L 

BRUSSELS 

0806 

C 

MAUI  ISLAND 

I 

43848 

J 

BUENOS  AIRES 

8817 

J 

OAHU- HONOLULU 

u 

43848 

K 

BUENOS  AIRES 

I 

9087 

J 

MIDWAY 

I 

43008 

K 

BULAWAYO 

965  1 

J 

TAIPEI 

1 

21523 

V 

BURG 

9667 

J 

KAO ENA 

38265 

A 

CAIRNS 

9724 

L 

HONG  KONG 

10591 

M 

CAIRO 

I0028 

J 

CALCUTTA 

10028 

J 

CALCUTTA 

10052 

J 

BANARAS 

836 

K 

CALI 

10060 

J 

LUCKNOW 

I 

40111 

J 

CAMPOS 

10132 

J 

AHMF.DABAD 

45459 

J 

CANBERRA 

10143 

J 

UDAIPUR 

I 

43858 

J 

CANUELAS 

10165 

J 

JAIPUR 

u 

3962 

J 

CAPE  VERDE  IS 

10177 

J 

AGRA 

I 

46738 

A 

CAPETOWN 

10187 

K 

NEW  DELHI 

u 

46738 

K 

CAPETOWN  K 

105  1 1 

K 

WADI  HALE A 

4306 

K 

CARACAS 

I 

10542 

K 

ASWAN 

I 

36479 

J 

CARAVEL AS 

I 

10552 

K 

LUXOR 

47503 

K 

CARMEN  DE  PATAGONES 

1059  1 

M 

CAIRO 

1 

32977 

J 

CAROLINA 

10909 

J 

AG AD  I R 

10937 

J 

CASABLANCA 

10937 

J 

CASABLANCA 

u 

17930 

N 

CASTIGLIONCEL 

10955 

J 

TANGIER 

u 

1 4375 

B 

CATANIA  B 

u 

10966 

K 

ROIA 

I 

14395 

N 

CETRARO 

10989 

K 

LISBON 

18070 

K 

CHATEAU  RENAULT 

11187 

J 

AZORES 

48732 

K 

CHRISTCHURCH 

1 1524 

J 

BERMUDA 

u 

38726 

J 

COCOS  1SL.  J 

u 

1 1524 

K 

BERMUDA 

1796  1 

J 

COLLE  ISARCO 

I 

1 3080 

A 

TOHOKU 

2969 

J 

COLOMBO 

131  10 

A 

KAGOSH I MA 

47557 

K 

COMODORO  RIVADAVIA 

13120 

A 

KUMAMOTO 

I 

293 

J 

CONAKRY 

u 

13130 

A 

KYUSHU 

1 

21552 

K 

COPENHAGEN 

13145 

J 

ITAMI 

u 

21552 

C 

COPENHAGEN  C 

13155 

C 

KYOTO 

439  1 4 

K 

CORDOBA 

I 

13159 

c 

TOKYO 

35769 

K 

DAR  ES  SALAAM 

CODE: 

I = > 

1 IGSN71  net  only,  0 *> 

(JAU  net  only,  NO 

CODE  = > common  to  both 
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TABLE  1 (continued) 


WORLD  NETS 

exc 

ludlng  NORTH  AMERICAN 

NETS  ( c o n t i nue  d ) 

Lis  t inn 

i n 

1GB  NUMBER  order 

List! ng  1 n 

ALPHABETICAL  order 

CODE 

IGB 

IGB 

U 

13159 

N 

TOKYO  N 

I 

38320 

A 

DARWIN 

U 

13276 

J 

SEOUL 

U 

38320 

J 

DARWIN  J 

1 

13707 

J 

MOHAN 

13708 

A 

DEHRA  DUN 

13700 

A 

DEHRA  DUN 

U 

14395 

L 

D 1 AMANTE 

1 

137  14 

J 

AMRITSAR 

3649 

J 

DOUALA 

I 

1 3049 

J 

KABUL 

48750 

D 

DUNEDIN 

1 

1395  1 

J 

TEHERAN 

I 

18153 

J 

EDINBURGH 

14  l 12 

K 

PORT  SAID 

u 

18153 

0 

ED  I N BURGH 

14135 

J 

BE  1 PUT 

21550 

P 

El  BY 

( 

14192 

J 

ANKARA 

3302 

J 

ENTEBBE 

1 4323 

A 

TR 1 POL  1 

I 

14374 

T 

ETNA  KM  15-16 

U 

14374 

P 

ETNA  1’ 

U 

14374 

P 

ETNA  P 

I 

14374 

T 

ETNA  KM  15-16 

14386 

J 

FALERNA  MARINA 

u 

14375 

B 

CATANIA  B 

25  1 75 

J 

FAUSKE 

I 

14375 

X 

S.  BERNARDO 

I 

17941 

F 

FERRARA 

u 

14385 

J 

GALATI  J 

I 

40178 

J 

FLORI ANOPOLIS 

I 

14385 

T 

Al.l  TERM}', 

25142 

R 

FORMOFOSS 

14386 

J 

EALERNA  MARINA 

32838 

J 

FORTALEZA 

u 

14395 

L 

D I AMANTE 

I 

21608 

0 

FRANKFURT 

I 

14395 

N 

CETRARO 

u 

21608 

P 

FRANKFURT  P 

14396 

J 

S.  LUC  IDO 

u 

21609 

T 

FULDA 

I 

14463 

J 

ALG I ERS 

u 

14385 

J 

GALATI  J 

14492 

J 

MALLORCA 

668 

J 

GEORGETOWN 

14503 

M 

MADRID 

I 

18154 

P 

GLASGOW 

16601 

J 

M1SAWA 

I 

36569 

J 

GO  1 ANA 

16631 

K 

SAPPORO 

I 

41792 

K 

CRAFTON 

1 

1665  1 

A 

WAKKANAI 

7485 

J 

GRAND  CANARY 

u 

17904 

J 

ANGRI  J 

5834 

N 

GUAM 

I 

17904 

P 

LI COLA 

I 

4495 

J 

GUANTANAMO 

u 

17905 

J 

PONTE  FARAONE 

33229 

K 

GUAYAQUIL 

I 

17905 

L 

BIVIO  GIUNGANO 

59520 

J 

HALLETT 

I 

179  12 

A 

ROME 

25101 

K 

HAMAR 

u 

179  12 

N 

ROME  N 

21639 

B 

HAMBURG 

17913 

N 

MINTURNO 

28603 

A 

HAMMERFEST 

17921 

J 

PODERE  SP I NET A 

D 

21629 

A 

HANNOVER  A 

I 

17930 

J 

OUERCETA 

I 

21629 

K 

HANOVER 

u 

17930 

N 

CAST  I CL  I ONCEL 

45196 

J 

HASTINGS 

1 

17940 

J 

LUZZARA 

U 

5295 

J 

HAWAII  ISLAND 

u 

17940 

P 

RICO 

21521 

J 

HELMSTEDT 

I 

17941 

E 

FERRARA 

U 

21562 

J 

1 ELS  I NGBORG  J 

17950 

J 

PERI 

I 

21562 

T 

HELS I NGOR 

17951 

G 

ROVERETO 

I 

25004 

A 

HELSINKI 

17961 

J 

COLLE  ISARCO 

u 

25004 

S 

HELSINKI  S 

u 

17971 

K 

INNSBRUCK 

I 

25229 

u 

HJERK1NN 

I 

17971 

X 

STAFFLACH 

V 

21581 

J 

HOGS  TORI’  3 

17972 

L 

NIEDERAUDORF 

9724 

L 

HONG  KONG 

u 

17981 

C 

MUNICH  C 

I 

6578 

J 

HYDERABAD 

I 

17981 

,1 

MUNICH 

u 

17971 

K 

INNSBRUCK 

17990 

B 

BAMBERC  SUD 

4040O 

K 

IQU1QUE 

u 

1799  1 

D 

NURNBERG 

I 

33233 

J 

IQUITOS 

I 

1799  1 

P 

NEUSES 

13145 

J 

IT  AMI 

18012 

J 

BARCELONA 

25087 

J 

IVALO 

18022 

.1 

PERPIGNAN 

10165 

J 

JAIPUR 

18030 

L 

TARBES 

43068 

L 

JOHANNESBURG 

1803  1 

J 

TOULOUSE 

I 

13849 

J 

KABUL 

I 

18033 

J 

NARBONNE 

9667 

J 

KADENA 

1 

1 8035 

C 

MARSEILLES 

131  10 

A 

KAGOSHIMA 

u 

18037 

J 

NICE 

I 

7228 

J 

KANO 

18040 

J 

AGEN 

u 

21619 

R 

KASSEL  OST 

18049 

R 

PI  ASTRA 

I 

45312 

J 

KEMPSEY 

18050 

J 

MONTIGNAC 

6952 

K 

KHARTOUM 

18059 

J 

MILAN 

43084 

K 

KIMBERLEY 

18060 

K 

POITIERS 

I 

4476 

J 

KINGSTON 

18070 

K 

CHATEAU  RENAULT 

u 

4476 

L 

KINGSTON  L 

ICSN7 1 net  only,  U =>  UAU  net  only,  NO  CODE  =>  common  to  both  nets 
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CODE:  I = > 


TABLE  I (continued* 


WORLD  NETS  excluding  NORTH  AMERICAN  NETS  (continued) 


L Is  t 1 ng 

In 

ICI1  NUMBER  order 

List in« 

1 n 

ALPHABETICAL  order 

I OB 

IGB 

CODE 

1 8082 

0 

PARIS 

35945 

M 

K 1 NS U ASA/ LEOPOLD V I LL 

I 

181  10 

A 

TEDDING  'ON 

263  1 

J 

KUALA  LUMPUR 

U 

18110 

J 

TEDDINCTON  J 

13120 

A 

KUMAMOTO 

I 

18153 

J 

EDINBURGH 

u 

2087 

J 

KWAJALE 1 N 

u 

18153 

0 

EDINBURGH 

13  155 

c 

KYOTO 

I 

18154 

P 

GLASGOW 

u 

13130 

A 

KYUSHU 

I 

18165 

J 

OBAN 

I 

36768 

A 

LA  BY/ 

21500 

"L 

BAD  NEUSTADT 

u 

36768 

J 

LA  PAZ  J 

I 

2 15  10 

A 

BAD  HARZRURG 

3663 

J 

LAGOS 

u 

215  10 

c 

BAD  HARZBURG 

46630 

J 

LA ! NGHURG 

21520 

C 

BRAUNSCHVE 1 < 

25165 

K 

LE I RJ  ORDF  ALL 

21521 

J 

HELMS' TEDT 

3609 

H 

LIBREVILLE 

I 

21528 

A 

POTSDAM 

I 

17904 

P 

LI  COLA 

I 

21523 

V 

BURG 

25  110 

P 

LILLEHAMMER 

u 

21523 

w 

BRANDENBURG 

36827 

K 

1.1MV 

21330 

L 

ST00KELD0 RF- F ACRE N BU 

109  89 

K 

LISBON 

I 

21540 

.1 

RIGK1.INC 

u 

43055 

B 

LOB ATS  1 

21550 

P 

F.  I BY 

u 

42952 

J 

LOURENCO  MARQ 

2 155  1 

J 

R 1 NESTED 

35983 

J 

LUANDA 

u 

21552 

C 

COPENHAGEN  C 

10060 

J 

LUCKNOW 

I 

21552 

K 

COPENHAGEN 

39458 

J 

LUSAKA 

u 

21362 

J 

HELS INGBORG  J 

1 

10532 

K 

LUXOR 

l 

21562 

T 

HELS1NC0R 

1 

17940 

J 

LUZZARA 

21563 

J 

VEINGE  KE. 

4 1819 

J 

MACKAY 

21571 

J 

S.  KRISTINA 

I 

64  30 

J 

MADRAS 

u 

2 1572 

J 

APELVIKSAAS  J 

14503 

M 

MADRID 

l 

21572 

K 

APELVIKSAAS 

25  1 3 1 

K 

MAERE 

u 

21581 

J 

HOGS TO RP  J 

25  1 5 3 

K 

MAJAVATN 

I 

21581 

U 

TANUM 

u 

2622 

J 

MALACCA 

21590 

K 

OSLO 

14492 

J 

MALLORCA 

21591 

J 

SVINESUNDE 

I 

33039 

J 

MANAUS 

I 

21397 

K 

STOCKHOLM 

6050 

L 

MAN  I LA 

u 

21604 

L 

BRUSSELS 

I 

59637 

J 

MARBLE  POINT  J 

I 

21604 

S 

BRUSSELS 

I 

18035 

C 

MARSEILLES 

1 

21608 

0 

FRANK!  URT 

I 

4 1752 

A 

MARYBOROUGH 

u 

21608 

P 

FRANKFURT  P 

u 

793 

K 

MA  TURIN 

I 

21609 

J 

BAD  HERSFELD 

8806 

C 

MAUI  ISLAND 

u 

21609 

T 

FULDA 

42707 

J 

MAURITIUS  ISLAND 

u 

216  19 

R 

KASSEL  OST 

u 

4296  1 

B 

MBABANE 

I 

216  19 

V 

MELS UNGEN- BEUERN 

I 

35783 

K 

MBEYA 

u 

21629 

A 

HANNOVER  A 

3846 

J 

MBOUR- DAKAR 

I 

21629 

K 

HANOVER 

59676 

C 

MCMURDO  SOUND 

I 

21629 

R 

SOLTAU 

865 

K 

MKDELL I N 

I 

21638 

K 

BREMEN 

45474 

M 

MELBOURNE 

21639 

B 

HAMBURG 

I 

216  19 

V 

MELSUNCEN-BEUERN 

21649 

F 

RENDSBURG 

21659 

J 

M I DDELF  ART 

21659 

J 

M1DDELF ART 

I 

9087 

J 

M 1 DVAY 

I 

21716 

P 

TVERAA 

18059 

.1 

MILAN 

21941 

K 

REYKJAVIK 

1 79  1 3 

N 

M IN TURN 0 

u 

22270 

.1 

SONDRF.STROMF  J 

16601 

J 

M 1 SAW A 

I 

25004 

A 

HELSINKI 

25  164 

K 

MO-  I- RAN  A 

u 

25004 

s 

HELSINKI  S 

I 

1 3707 

J 

MOHAN 

I 

25043 

J 

OULU 

u 

35749 

D 

MOMBASA 

25065 

J 

ROVANEIMI 

1 

260 

K 

MONROV I A 

25087 

.1 

IVALO 

1 

4487 

,| 

MONTEGO  HAY 

I 

25090 

J 

SORKJOSEN 

43846 

K 

MONTEVIDEO 

25093 

J 

ALTA 

18050 

J 

MONTICNAC 

25101 

K 

HAMAR 

I 

35737 

K 

MOSR  1 

25  1 10 

P 

LILI.EHAMMER 

I 

4 1909 

J 

MT.  ISA 

25  120 

J 

SOKNEDAL 

1 

17981 

J 

MUN ICH 

25130 

L 

TRONDHEIM 

u 

17981 

C 

MUNICH  C 

23  131 

K 

MAERE 

I 

35716 

A 

NAIROBI 

25142 

R 

FORMOFOSS 

II 

35716 

N 

NAIROBI  N 

25  1 43 

J 

VE I SKILLS 

37977 

J 

NANDI-EIJ I ISLAND 

CODE: 

: I => 

IGSN71  net  only,  U =>  UAU 

net  only,  NO 

CODE  =>  common  to  both  i 
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TABLE  1 (continued) 


WORLD  RETS  excluding  NORTH  AMERICAN  NETS  (continued) 


Listing  In  IGB  NUMBER  order 
IGB 

CODE  NUMBER  NAME 

25153  K MAJAVATN 
I 25163  J SKAMDAL 
25164  K MO- I -RAN A 
25  165  K LE1RJORDFALL 

25174  .1  BODO 

25175  J FAUSKE 
251B7  K NARVIK 

25198  K BARDUEOSS 

25199  J TROMSO 
25219  Q VINSTRA 

U 25229  L Ol’PDAL  L 
I 25229  U HJERKINN 
25968  K THULE 
28603  A I1AMMERFEST 
I 29522  J ALERT 

32674  J ASCENSION  ISLAND 
32838  J FORTALEZA 
U 32084  J RECIFE  J 
1 32884  L RECIFE 

32918  L BELEM 
I 32977  J CAROLINA 
I 33039  J MANAUS 
1 33134  J TEFE 

33208  K GUI  TO 
33229  K GUAYAQUIL 
1 33233  J 1QU1T0S 

33341  K TALARA 
1 35716  A NAIROBI 

IT  35716  N NAIROBI  N 
I 35737  K MOSHI 
U 35749  D MOMBASA 

35769  K DAR  ES  SALAAM 
I 35783  K MBEYA 

35945  M K1NSHASA/LE0P0LDVILL 
35983  J LUANDA 
I 36428  J SALVADOR 
I 36479  J CAKAVELAS 
I 36508  J PORTO  NATIONAL 
I 36557  J BRAZIL! A 
I 36569  J GO  I ANA 
I 36593  J BELO  HORIZONTE 
I 36768  A LA  PAZ 
U 36768  .1  LA  PAZ  J 
U 36773  J SANTA  CRUZ 
36827  K LIMA 
36861  K AREGUIPA 
36880  K ARICA 
U 37579  B TAHITI 
U 37841  J PAGO  PAGO 
I 37841  L PAGO  PAGO 

37977  J NANDI -FIJI  ISLAND 
38265  A CAIRNS 
38296  N TOWNSVILLE 
I 38320  A DARWIN 
U 38320  J DARWIN  J 
U 38726  J COCOS  ISL.  J 
39371  M SALISBURY 
39428  K NDOLA 
39458  J LUSAKA 
39475  K VICTORIA  FALLS 
I 39525  J NOVA  LISBOA  J 
I 39543  J SA  DA  BANDEIRA 
40100  .J  VITORIA 
I 40 1 I l J CAMPOS 


Listing  In  ALPHABETICAL  order 
IGB 

CODE  NUMBER  NAME 

I 18033  J NARBONNE 
25187  K NARVIK 
39428  K NDOLA 
I 17991  P NEUSES 
10187  K NEW  DELHI 
I 7232  J NIAMEY 

U 18037  J NICE 

17972  L NIEDERAUDORF 

3885  J NOUAKCHOTT 

1 39525  J NOVA  LISBOA  J 

U 17991  D NURNBERG 

8817  J OAHU- HONOLULU 
I 18165  J OBAN 
U 25229  L OP P DAL  L 
40334  K ORAN 
21590  K OSLO 
I 25045  J OULU 
I 37841  L PAGO  PAGO 
U 37841  J PAGO  PAGO 
655  J PARAMARIBO 
18082  0 PARIS 
I 43812  J PELOTAS 
2650  J PENANG 
17950  J PERI 
18022  J PERPIGNAN 
I 45715  P PERTH 
U 45715  A PERTH  A 
18049  R PI  ASTRA 
I 43039  K P I ETERSBURG 
1792 1 J PODERE  SPINETA 
18060  K POITIERS 
U 17905  J PONTE  FARAONE 
826  K POPAYAN 
I 4482  J PORT  AU  PRINCE 
7407  J PORT  ETIENNE 

4301  J PORT  OF  SPAIN 

14112  K PORT  SAID 

6997  K PORT  SUDAN 

I 43801  J PORTO  ALECRE 
I 36508  J PORTO  NATIONAL 
I 21523  A POTSDAM 
43058  A PRETORIA 
47575  K PUERTO  DESEADO 
I 47612  J PUERTO  MONTT 

51108  K PUERTO  SANTA  CRUZ 
51230  L PUNTA  ARENAS 
I 17930  J QUERCETA 
33208  K QUITO 
4387  K RAMEY 
6366  K RANGOON 
I 32884  L RECIFE 
U 32884  J RECIFE  J 
21649  F RENDSBURG 
21941  K REYKJAVIK 
I 21540  J RICKLING 
U 17940  P RICO 

21551  J RINGSTED 
43934  K RIO  CUARTO 
40123  A RIO  DE  JANEIRO 
51119  K RIO  GALLEGOS 
51137  L RIO  GRANDE 
41730  K ROCKHAMPTON 
I 17912  A ROME 
U 17912  N ROMF.  N 


CODE 


TABLE  1 (continued) 


WORLD  RETS  excluding  NORTH  AMERICAN  NETS  (continued) 


List  Inp; 

in 

ICB  NUMBER  order 

L 1 h 1 1 ng 

i n 

ALPHABETICAL  order 

ICB 

IGB 

CODE 

CODE 

40123 

A 

RIO  DE  JANEIRO 

43920 

K 

ROSARIO 

40136 

J 

SAO  PAULO 

U 

10966 

K 

ROTA 

I 

40178 

J 

FLORIANOPOLIS 

25065 

J 

ROVANE I MI 

40257 

J 

ASUNCION 

1795  1 

G 

ROVERETO 

40334 

K 

ORAN 

I 

14375 

X 

S.  BERNARDO 

40345 

K 

SALTA 

21571 

J 

S.  KRISTINA 

40365 

L 

TUCUMAN 

14396 

J 

S . LUC  I DO 

40374 

K 

SANTIAGO  ESTERO 

I 

39543 

J 

SA  DA  BANDEIRA 

40400 

K 

IQUIOUE 

6206 

J 

SAIGON 

40420 

K 

TOCOFILLA 

3937  1 

M 

SALISBURY 

40430 

K 

ANTOFAGASTA 

40345 

K 

SALTA 

41730 

K 

ROCKHAMPTON 

I 

36428 

J 

SALVADOR 

I 

41752 

A 

MARYBOROUGH 

4386 

J 

SAN  JUAN 

U 

41752 

J 

MARYBOROUGH 

47597 

K 

SAN  JULIAN 

41773 

J 

BRISDANE 

U 

36773 

J 

SANTA  CHUZ 

I 

41792 

K 

GRAFTON 

I 

44030 

A 

SANTIAGO 

41819 

J 

MACKAY 

40374 

K 

SANTIAGO  ESTERO 

1 

41909 

J 

MT.  ISA 

u 

44030 

K 

SANTIAGO  K 

I 

41933 

J 

ALICE  SPRINGS 

40136 

J 

SAO  PAULO 

42707 

J 

MAURITIUS  ISLAND 

16631 

K 

SAPPORO 

U 

42952 

J 

LOURENCO  M.ARQ 

u 

13276 

J 

SEOUL 

U 

4296  1 

B 

MBABANE 

26  13 

A 

SINGAPORE 

I 

43008 

K 

BULAWAYO 

I 

25  163 

J 

SKAMDAI. 

I 

43039 

K 

P 1 ETERSBURG 

25120 

J 

SOKNEDAL 

u 

43055 

B 

LOBATS I 

I 

21629 

R 

SOLTAU 

43058 

A 

PRETORIA 

u 

22270 

J 

SONDRESTROMFJ 

43068 

L 

JOHANNESBURG 

2670 

J 

SONGKHLA 

43084 

K 

KIMBERLEY 

F 

25090 

J 

SORKJOSF.N 

I 

43B01 

J 

PORTO  ALEGRE 

I 

4374 

J 

ST.  CROIX 

I 

43812 

J 

PELOTAS 

1 

4341 

J 

ST.  LUCIA 

43846 

K 

MONTEVIDEO 

I 

17971 

X 

STAFFLACH 

I 

43848 

J 

BUENOS  AIRES 

I 

21597 

K 

STOCKHOLM 

u 

43848 

K 

BUENOS  AIRES 

21530 

L 

STOOKELDORF-FACKENBU 

I 

43858 

J 

CANUELAS 

2159  1 

J 

SV 1 NESUNDE 

439  1 4 

K 

CORDOBA 

45331 

J 

SYDNEY 

43920 

K 

ROSARIO 

u 

37579 

B 

TAHITI 

43934 

K 

RIO  CIJARTO 

965  1 

J 

TA I PE  I 

43982 

K 

BAHIA  BLANCA 

33341 

K 

TALARA 

I 

44030 

A 

SANTIAGO 

10955 

J 

TANG  I ER 

u 

44030 

K 

SANTIAGO  K 

I 

21581 

Q 

TANUM 

u 

44031 

K 

VALPARAISO 

18030 

L 

TARBES 

u 

45  164 

0 

AUCKLAND 

I 

181  10 

A 

TEDDI NGTON 

45196 

.1 

HASTINGS 

u 

181  10 

J 

TEDDI NGTON  J 

I 

45312 

J 

KEMPSEY 

I 

33134 

J 

TEFE 

45331 

J 

SYDNEY 

1 

1395  1 

J 

TEHERAN 

45459 

J 

CANBERRA 

I 

6956 

J 

TESSENE I 

I 

45466 

K 

ALBURY 

25968 

K 

THULE 

45474 

M 

MELBOURNE 

40420 

K 

TOGOP I LLA 

u 

45715 

A 

PERTH  A 

I 

13080 

A 

TOHOKU 

I 

45715 

P 

PERTH 

1 

13159 

C 

TOKYO 

46603 

J 

BRISTOWN 

u 

13159 

N 

TOKYO  N 

46622 

J 

BEAUFORT  WEST 

1 803  1 

J 

TOULOUSE 

46630 

J 

LAINGBURG 

38296 

N 

TOWNSV I LLE 

1 

46738 

A 

CAPETOWN 

47535 

K 

TRELEW 

u 

46738 

K 

CAPETOWN  K 

1 4323 

A 

TRIPOLI 

47503 

K 

GARMF.N  DE  PATAGONES 

25  199 

J 

TROMSO 

47535 

K 

TRELEW 

25130 

L 

TRONDHEIM 

47557 

K 

COMODORO  RIVADAVIA 

40365 

L 

TUCUMAN 

47575 

K 

PUERTO  DESEADO 

I 

21716 

P 

TVERAA 

47597 

K 

SAN  JULIAN 

10143 

J 

UDAIPUR 

1 

476  12 

J 

PUERTO  MONTT 

51148 

1, 

USHUA 1 A 

48714 

K 

WELL  I NGTON 

u 

4403  1 

K 

VAL.PARAISO 

48732 

K 

CHRISTCHURCH 

21563 

J 

VEINGE  KE. 

48750 

D 

DUN  El)  1 N 

25  1 43 

J 

VEISKILLE 

: I =>  ICSN7I  net  only,  0 - > UAU  net  only,  NO  CODE  =>  common  to  both  nets 

2 1 


CODE 


TABLE 


WORLD  NETS  excluding;  NORTH  AMERICAN 


Listing;  In  IGB  NUMBER  order 
IGB 

CODE  NUMBER  NAME 

51108  K PUERTO  SANTA  CRUZ 
51119  K RIO  GALLEGOS 
5 1 137  L RIO  GRANDE 
51148  L USHUAIA 
51230  L PUNTA  ARENAS 
59520  J HALLETT 
I 59637  J MARBLE  POINT  J 
59676  C MCMURDO  SOUND 


1 (contlnned) 


NETS  (continued) 


Listing;  In  ALPHABETICAL  order 
IGB 

CODE  NUMBER  NAME 

39475  K VICTORIA  FALLS 
25219  O VINSTRA 
40100  J VITORIA 
10511  K WADI  HALFA 
I 5696  N WAKE  ISLAND 
U 5696  J WAKE  ISLAND  J 
I 16651  A WAKKANAI 
48714  K WELLINGTON 


CODE:  I =>  IGSN71  net  only,  U =>  UAU  net  only,  NO  CODE  =>  common  to  both  nets 


f 


rr 


4.  Bogota 

5.  Bodo 

6.  Maryborough 


884  K 
25  174  J 
41752  A. 

During  and  after  these  analyses  were  done,  new  absolute  measure- 
ments were  made  with  Italian  absolute  apparatus.  All  new  absolute 
measurements  could  not  be  included  in  the  new  analyses 
because  they  were  not  made  at  1GSN  7 1 stations  and  we  did  not 
have  information  about  ties.  New  absolute  measurements  were  con- 
sidered at  the  following  stations: 


Teddington 

18  1 10 

A 

Rome 

179  12 

A 

Hammerfest 

28603 

A 

Helsinki 

25004 

A 

M unic  h 

1798  1 

A 

Copenhagen 

2 1552 

K 

S.  Bernardo 

14375 

X 

B run  sc  hweig 

2 1520 

C 

Hamburg 

2 1639 

B 

Pa  ri  s 

18082 

O. 

We  did  not  have  good 

info  rmation 

measurements,  but  we  used  0.02  mgal.  After  including  these 
measurements  in  IGSN  7 1 net,  the  stations  where  new  absolute 
measurements  should  be  made  came  in  the  following  order: 


1. 

Nai robi 

357  16 

A 

2. 

New  Delhi 

10187 

K 

3. 

Rio  Gallegos 

5 1119 

K 

4. 

Bogota 

884 

K 

5. 

Rockha  mp  ton 

4 1730 

K 

6. 

McMurdo  Sound 

59676 

C 

In 

the  adjustment  of 

IGSN  7 1 

on' 

solved  for  the  scales  of  the  gravimeters.  We  had  in  hand  a 


24 


1 


variance-covariance  matrix  for  the  gravity  values  solved  under  the 
AFCRL  Contract  No.  F 19628 -68 -C -0335  (UAU-net),  which  included 
also  the  second  order  corrections  to  scales  of  the  gravimeters. 

We  did  the  analyses  for  this  net  which  included  372  gravity  stations 
distributed  around  the  world,  listed  in  Table  1 and  shown  in  Figure  2. 

Two  separate  analyses  were  done:  1)  original  net  and  2) 
original  net  plus  ten  new  absolute  measurements  li  sted  above.  In  the 
order  of  preference  six  new  sites  for  future  absolute  measurements 
were  selected  as  follows: 

1)  Original  net 


1. 

Panama 

889  A 

2. 

Thule 

25968  K 

3. 

Nai  robi 

35176  N 

4. 

Wa  shington 

11687  M 

5. 

Sydn  ey 

45331  J 

6. 

Bodo 

25174  J. 

With  new  absolute  measurements 

1. 

Paso  De  Cortes 

4698  A 

2. 

Singapore 

2613  A 

3. 

Thule 

25968  K 

4. 

Azores 

11187  J 

5. 

Nairobi 

357  16  N 

6. 

Buenos  Aires 

43848  K. 

The 

se  two  nets;  IGSN 

7 1 and  the  UAU-net  do  not  include  the 

same  observations  and  not  necessarily  the  same  statiors  as  seen  in 
Table  1.  The  latter  one  included  the  same  pendulum  and  absolute 
measurements  as  IGSN  7 1 but  from  gravimeter  ties  only  those 
which  were  made  with  La  Coste-Romberg  gravimeters  were  included 
in  the  adjustment  of  the  UAU-net.  When  we  looked  at  the  distribu- 
tion of  the  selected  stations  for  new  absolute  measurements,  it  is 
our  opinion  that  this  latter  selection  might  be  a better  one.  It  also 
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is  supposed  to  control  better  the  second  order  influence  in  the 
scale  factors. 

More  detailed  analyses  showed  that  after  the  new  absolute 
measurements  in  Europe,  new  proposed  absolute  measurements 
contributed  hardly  anything  to  the  European  net.  The  contributions 
of  new  absolute  measurements  have  become  more  continental  and 
local  rather  than  global.  It  lias  become  obvious  that  we  do  not 
gain  much  through  old  gravimeter  connections  between  continents. 

If  we  wish  to  improve  world  wide  net,  the  best  improvements  are 
coming  from  new  absolute  measurements  on  different  continents. 

It  also  shows  that  the  best  approac  h to  improve  IGSN  7 1 net  is  to  do 
readjustments  of  continental  or  local  nets  rather  than  a new  adjust- 
ment of  the  global  net.  Even  new  gravimeter  ties  between  continents 
with  the  currently  available  instruments  will  not  change  much  the 
situation. 

4.6.2  African  network 

Outside  of  the  contract  we  made  analyses  of  African  portion 
of  IGSN  7 1 net  as  a favor  to  African  nations,  who  were  planning 
to  establish  new  absolute  gravity  sites  in  Africa.  We  selected  the 
stations  from  IGSN  7 1 net  which  fell  in  Africa  plus  Paris  and  Rome, 
where  new  absolute  measurements  had  been  made  at  that  time 
(March,  1977).  We  wanted  to  see  if  these  same  selection  criteria 
could  be  used  for  a smaller  part  of  the  world  net.  The  preferred 
order  of  new  absolute  measurements  came  as  follows: 


1. 

Kin  sha  sa 

35945 

M 

2. 

Nairobi 

357  16 

A 

3. 

Mbour-Daka  r 

3846 

J 

4. 

Lusaka 

39458 

J 

6. 

Ca  sablanca 

10937 

J 

27 

6. 

Luxo  r 

10552 

K 

7. 

Douala 

3649 

J 

8. 

Nova  Lisboa 

39575 

J 

9. 

Beaufort  West 

46622 

J 

10. 

Da  r Es  Salaam 

35769 

K 

11. 

A sma  ra 

6958 

K 

12. 

Port  Etienne 

7407 

J 

13. 

Damako 

3728 

J 

14. 

Bulawayo 

43008 

K 

15. 

Monrovia 

260 

K 

16. 

Cai  ro 

1059  1 

M. 

This  net  included  only  the  stations  which  were  in  IGSN  7 1 net. 
Therefore,  it  might  not  be  the  best  for  Africa  as  far  as  a new 
base  station  net  is  concerned.  It  might  be  necessary  to  establish 
more  dense  net  taking  into  consideration  existing  gravity  measure- 
ments and  a future  use  of  the  net.  However,  analyses  of  the 
results  indicated  that  the  program  was  making  the  logical  choices 
as  far  as  the  IGSN  7 1 improvements  were  concerned. 

4.6.3  North  American  and  U.S.  Gravity  Base  Station  nets 

During  the  summer  of  1977  the  negotiations  were  going  on  to  bring 
Italian  absolute  apparatus  to  the  United  States  for  inter-comparison 
purposes.  We  were  asked  by  Project  Monitor,  Bela  Szabo  to 
analyze  the  U.S.  portion  of  IGSN  7 1 net  for  the  preferred  locations 
of  new  absolute  measurements.  Three  versions  were  studied: 

1.  Stations  in  North  America. 

2.  Stations  in  North  America  but  selected  stations  for 
absolute  measurements  to  be  located  in  the  United 
States,  excluding  Alaska. 

3.  Stations  in  the  United  States  excluding  Alaska  and 
Hawaii. 
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All  of  the  above  three  alternatives  were  analyzed  in  two  ways: 

a)  considering  stations  in  IGSN  7 1 and  the  variance- 
covariance  matrix  from  the  solution,  where  only 
linear  correction  term  to  calibration  of  gravimeters 
were  included; 

b)  considering  the  UAU-network  and  the  variance- 
covariance  matrix  obtained  from  the  solution,  which 
included  linear  and  second  order  correction  terms  to 
the  calibration  of  gravimeters. 

In  Table  1 the  stations  belonging  to  the  U.S.  networks  are  identi- 
fied as  well  as  other  stations  belonging  to  the  North  American  nets. 

All  of  the  above  mentioned  alternatives  were  examined. 

During  the  analyses  not  only  the  preferred  order  was  determined 
for  the  new  absolute  measurements,  but  several  quantities  were 
computed,  such  as  the  trace  of  variance-covariance  matrix,  the 
partial  trace  of  the  same  matrix  as  explained  earlier,  new 
variances  of  the  gravity  values  of  the  stations,  average  variances 
for  the  stations,  changes  in  variances  and  percentage  changes  of 
variances.  These  quantities  were  computed  after  each  cycle  for 
each  station  included  in  the  solutions.  In  the  following  three  tables, 

2,  3 and  4,  the  solutions  using  122  stations  in  North  America  are 
summarized.  The  selected  stations  are  given  in  order  of  preference. 
In  these  analyses  it  was  assumed  that  the  new  absolute  measurements 
have  an  accuracy,  a = 20  pgal.  Table  2 gives  results  for  both  of 
the  networks  and  the  average  variances  of  stations  in  North 
America  at  the  beginning  and  after  each  addition  of  the  absolute 
measurement.  In  these  solutions  there  were  no  preselected  stations 
for  absolute  measurements.  They  were  "free  solutions." 

Table  3 gives  the  selections  of  preferred  stations  considering 
the  effect  of  new  absolute  measurements  in  the  whole  North  American 
nets,  but  constraining  the  selected  stations  to  be  located  inside 
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Table  2 

North  American  Base  Station  Nets 
F ree  Solution 


Accuracy  of  Absolute  Measurements  a - 20  jjgal 


122 

IGSN  7 1-Net 

, Stations,  a 2 = 711  ^gal2 

' 

122 

UAU-Net 

1 Stations,  5 2 = 1776  jg al2 

Selected  Stations 

mm 

Selected  Stations 

ct  2 Mgal2 

1. 

Mexico  City 

626 

1. 

Monterrey 

1255 

2. 

Hall  Beach 

588 

2. 

Point  Barrow 

808 

3. 

Dalla  s 

1 562 

3. 

Paso  De  Corte 

684 

4. 

Jacksonville 

545 

4- 

Wa  shington 

621 

5. 

Edmonton 

i 531 

! 5. 

Great  Falls 

591 

6. 

San  Jose 

1 5 18^  ^ m 

i‘^6. 

Resolute  Bay 

j 56 1 

7. 

Mould  Bay 

; 488 

1 

7. 

I 

San-M^ranci  sco 

544 

; 

8. 

Minneapoli  s 

479 

! 

8. 

D env  e r 

! 

534 

1 

Table  3 

North  American  Base  Station  Nets 
Free  Solutions  but  Stations  to  be 
Selected  in  U.  S.A. 


Accuracy  of  Absolute  Measurements  a = 20  jxgal 


IGSN  7 1-Net 

122  Stations,  5- 2 = 711  ;jgal2 

UAU-Net 

122  Stations,  52  = 1776  /ugal2 

Selected  Stations 

mm 

Selected  Stations 

a2  jigal2 

1. 

Miami 

639 

1. 

Miami 

1274 

2. 

San  Antonio 

609 

2. 

Denver 

1187 

3. 

Great  Falls 

588 

3. 

Sea  ttle 

1122 

4. 

Orlando 

574 

4. 

Miami 

1078 

5. 

Minneapoli  s 

562 

5. 

Wa  shington 

1044 

6. 

Albuquerque 

552 

6. 

El  Pa  so 

102  1 

7. 

Seattle 

543 

7. 

Caribou 

999 

8. 

Loui  svill  e 

534 

8. 

San  F ranci  sco 

987 

30 
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Table  4 

North  American  Base  Station  Nets 
Denver  and  Boston  Preselected  Stations 
All  Stations  Selected  in  the  U.S.A. 


Accuracy  of  Absolute  Measurements  a - 20  ^gal 


IGSN  7 1-Net 
122  Stations,  a2  = 711 

Ugal2 

12 

UAU-Net 

2 Stations,  ~2  - 1776  \jg al2 

Selected  Stations  i 

. 

~2  Mgal2 

Selected  Stations 

[ ftgal2 

r 

D en  v e r 

Denv  er 

| 

Bo  ston 

6 30 

Bo  ston 

1337 

1. 

Miami 

597 

1. 

Miami 

1143 

2. 

San  Antonio 

580 

2. 

Miami 

j 1097 

3. 

Seattle 

566 

3. 

1 

Sea  ttle 

1052 

4. 

Orlando 

556 

4. 

1 

C hey enne 

1026 

5. 

Bi  smarck 

547 

5. 

El  Pa  so 

i 1005 

6. 

El  Paso 

5 38 

6 . 

Seattle 

986 

Table  5 

U.  S.  Base  Station  Nets 
F ree  Solution 


Accuracy  of  Absolute  Measurement  a - 20  figal 


83 

IGSN7  1-Net 

Stations,  s2  = 576  /jgal2 

77 

UAU-Net 

Stations,  a2  - 10  30 

ggal2 

Selected  Stations 

Selected  Stations 

i 

5 2 /ugal2 

i. 

Miami 

499 

i. 

Houston 

631 

2. 

Da  11a  s 

466 

2. 

D env e r 

555 

3. 

San  Francisco  ; 

446  j 

3. 

Wa  shington 

5 19 

4. 

Charleston 

430 

4. 

Miami 

485 

5. 

Loui  sville 

413 

5. 

Great  Falls 

| 465 

6. 

Bi  smarck 

i 40 1 

I 

6 . 

O rlando 

! 452 

l 

7. 

Albuquerque 

392 

i 

7 • 

Madi  son 

440 

8. 

Orlando 

385 

1 8‘ 
i 

Albuquerque 

j 432 

of  the  United  States  excepting  Alaska.  Table  4 gives  the  solutions 
for  the  North  American  nets  with  same  constraints  as  in  Table  3 
and  in  the  addition  the  condition  that  Boston  and  Denver  have  been 
preselected  to  be  the  sites  of  the  new  absolute  measurements. 

Tables  5-9  give  results  for  various  situations  in  the  U.S. 
base  station  nets.  We  have  first  two  solutions  for  the  free  choice 
of  the  stations  without  any  preselected  stations,  given  in  Table  5 and 
6.  The  only  difference  between  these  solutions  is  that  in  Table  5 
the  accuracy  of  absolute  measurements  is  assumed  to  be  a = 20  /jga.1 
and  in  Table  6 -i  = 10  ugal.  The  orders  of  selections  are  different 
but  the  average  variances  have  not  improved  much  - only  13%  even 
though  the  accuracy  of  absolute  measurements  has  been  improved 
50%.  This  small  return  from  improvement  of  accuracies  of  abso- 
lute measurements  points  out  that  the  gravity  differences  between 
stations  must  be  measured  more  accurately  in  order  to  benefit  fully 
from  the  improved  accuracies  of  absolute  measurements.  This 
additional  accuracy  from  absolute  measurements  will  improve  only 
local  situations  or  those  stations  which  are  tied  more  accurately 
than  average  to  the  stations  where  absolute  measurements  are  made. 

It  is  interesting  to  note  that  if  new  absolute  measurements 
were  done  at  eight  sites,  then  there  would  not  be  much  difference 
in  IGSN  7 1,  if  all  stations  were  selected  freely  in  the  North  America 
or  in  the  U.S. A.  There  is  even  less  difference  if  the  sites  were 
freely  selected  in  the  U.S. A.  or  two  of  the  stations  were  preselected 
in  the  same  area.  However,  if  we  take  the  UAU-network  which 
includes  the  second  order  terms,  there  is  a large  difference  in  the 
average  variances  if  the  station  selection  is  limited  f>  the  area  of 
the  U.  S.A.  , but  not  much  effect  is  seen  by  preselecting  two  of  the 
stations  as  compared  with  the  "free  solution"  in  the  U.S.A.  We 
can  also  see  that  after  4-5  new  absolute  measurements  the  gain  in 
the  whole  U.S.  net  is  not  much  - the  improvements  will  be  more 
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of  a local  nature.  Tables  7-9  give  the  U.S.  networks  after  2,  4 
and  5 stations  have  been  preselected.  It  is  again  interesting  to 
see  that  a reasonable  preselection  of  stations  does  not  influence 
much  the  average  variances  in  the  whole  net.  In  Table  9 the 
average  variance  in  UAU-net  is  even  smaller  than  in  Table  5, 
which  is  a "free  solution."  We  have  to  remember  that  the 
selections  of  the  stations  have  been  made  using  "the  partial  trace" 
of  va  rianc e-cova rianc e matrices  in  order  to  minimize  the  influence 
of  local  stations;  therefore,  the  full  trace  of  the  variance-covariance 
matrix  is  not  necessarily  minimum  for  the  preferred  choices. 

It  is  clear  that  if  we  make  more  than  six  new  absolute 
measurements,  we  do  not  gain  much  as  far  as  the  current  base 
station  networks  are  concerned.  Larger  improvements  can  be 
expected  if  new,  more  accurate  measurements  of  gravity  differences 
between  the  stations  are  established. 

Since  these  studies  were  concluded  and  informally  reported 
to  Proj  ect  Monitor,  the  new  absolute  measurements  have  been 
carried  out  at  Bedford,  Denver,  Bismarck,  Miami,  San  Francisco 
and  Alamogordo,  which  corresponds  about  the  situation  given  in 
Table  8 after  two  selections.  Therefore,  it  is  appropriate  to  give 
expected  variances  to  all  stations  included  in  the  U.S. A.  portion 
of  the  UAU-net.  The  old  variances  and  the  new  ones  are  given  in 
Table  10.  As  we  can  see  in  the  UAU-net  there  are  five  stations 
which  are  poorly  tied  to  any  other  stations,  namely  Tampa,  Corpus 
Christy,  San  Diego,  Norton  AFB  and  Portland.  If  a full  trace 
would  have  been  used  in  selection,  these  stations  would  possibly 
have  come  up  as  the  first  choices,  but  the  whole  net  would  not 
have  been  improved  much.  The  variances  of  these  stations  are 
keeping  the  average  variance  large. 

For  easy  reference  the  changes  in  the  variance  are  given  in 
Table  11.  This  table  clearly  reflects  that  the  last  two  new 
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Table  6 

U.  S.  Base  Station  Nets 
F ree  Solution 


Accuracy  of  Absolute  Measurements  -r  = 10  ugal 


83 

1GSN  71-Net 

Stations,  d2  = 576  jugaP 

77 

UAU-Net 

Stations,  nl  - 1030  jigaP 

Selected  Stations 

c2  Mgal2 

Selected  Stations 

a1  jugal2 

1. 

Dalla  s 

450 

l 

Denver 

5 15 

2. 

Orlando 

419 

2. 

1 

Miami 

458 

3. 

Minneapoli  s 

396 

3. 

Great  Falls 

418 

4. 

Loui  sville 

378 

4 

Wa  shington 

405 

5. 

Charleston 

365 

5. 

San  Antonio 

396 

6. 

San  Francisco 

355 

1 

6. 

Bo  ston 

390 

7 o 

Minot 

344 

7. 

Kansas  City 

382 

8. 

Albuquerque 

336 

8. 

L 

O rlando 

377 

Table  7 

U.  S.  Base  Station  Nets 
Boston  and  Denver  Preselected  Stations 

Accuracy  of  Absolute  Measurements  r*  - 20  ygal 


IGSN  71-Net 

83  Stations,  a2  = 576  jugal2 

77 

UAU-Net 

Stations,  52  = 1030  jLgal2 

Selected  Stations 
! 

a 2 jugal2 

Selected  Stations 

a 2 j£al2 

Denver 

Denver 

Bo  ston 

486 

Bo  ston 

609 

i. 

Miami 

453 

1. 

Miami 

512 

1 

2. 

Albuquerque 

434 

2- 

Denver 

486 

3. 

C ha  rleston 

4 19 

3. 

O rlando 

469 

4. 

1 

Bi  smarck 

407 

4. 

Great  Falls 

452 

5. 

Loui  sville 

393 

5 

San  Antonio 

| 442 

1 

6. 

Dalla  s 

386 

6. 

Madi  son 

432 

! 
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Table  8 

U . S.  Base  Station  Nets 
Boston,  Denver,  Albuquerque,  and  Bismarck 

Preselected  St  tions 


Accuracy  of  Absolute  Measurements  - 20  ypal 


IGSN  71-Net 

8 3 Stations,  a2  = 576  yg al 2 

UAU-Net 

77  Stations,  ^2  1030  /jgal' 

Selected  Stations 

- qgal2  T 

Selected  Stations 

^gal“ 

Denver 

I ) i nv e r 

Bo  ston 

Bio  ston 

Bi  sma  rck 

Albuquerque 

Albuquerque 

! 440 

! 

Bi  sma  rck 

505 

1 

1.  Jacksonville 

420 

1.  Miami 

468 

2.  Louisville 

405 

2.  San  Francisco 

458 

3.  Miami  ! 

! 395 

\ ! 

3.  Orlando 

! 

44  1 

4.  Charleston 

t 

00 
— J 

1.  Madison 

i 432 

Table  9 

U.  S.  Base  Station  Nets 

Boston,  Denver, 

Albuquerque,  Bismarck 

and  Columbus 

Preselected  Stations 

Accuracy  of  Absolute  Measurements  a = 20  ygal 

IGSN  7 1-Net 

UAU-Net 

83  Stations,  J2  ~ 576  ^gal2 

77  Stations,  52  = 1030  ^gal2 

Selected  Stations 

1 o2  Ugal2 

Selected  Stations 

: 5 2 juga-12 

Denver 

[ 

Denver 

i 

• 

i 

Bo  ston 

i 

Bo  ston 

Albuquerque 

Albuquerque 

i 

! 

j 

Bi  sma  rck 

Columbus 

( 

j 

Columbus 

j 425 

Bi  sma  rck 

1 489 

| 

1.  Jacksonville 

407 

j 1.  Miami 

| 453 

2.  Miami 

j 397 

2.  Orlando 

i 440 

3.  Charleston 

388 

i 3.  San  Francisco 

I 

430 

L 


A 
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TABLE  10  - VARIANCES  IN  PIIRCOCAL  SGUARED 
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absolute  measurements  do  not  improve  the  whole  net  much  and  the 
improvements  are  more  local  than  at  the  beginning.  Another 
interpretation  could  be  that  the  correlation  is  getting  smaller  and 
accuracies  of  measured  gravity  differences  between  stations  are 
not  good  enough  to  transfer  information  from  new  absolute  measure- 
ments very  far  through  the  net.  If  we  leave  out  5 stations  in  the 
net,  which  have  the  variances  above  1000  ugal  then  the  average 
variance  after  two  selections  is  259  u,gal2  or  n - 16  pgal,  and  after 
four  selections  237  q,gal2  or  a = 15  p,gal.  The  gained  accuracy  from 
two  additional  absolute  measurements  would  not  be  significant. 

4.  7 Conclusion s 

The  original  scope  of  the  study  was  to  determine  the  pre- 
ferred locations  for  new  absolute  measurements  to  improve  IGSN  7 1 
global  net.  This  was  accomplished  in  the  early  part  of  the  contract 
as  reported  in  Chapter  6,  section  6.  1.  Six  stations  were  suggested 
as  the  preferred  sites  for  new  absolute  measurements.  It  was 
noted  that  the  last  two  new  measurements  made  contributions  more 
to  the  local  area  than  to  the  global  net.  During  the  work  new 

absolute  measurements  were  made  at  ten  more  sites  in  Europe.  A 
new  selection  of  the  preferred  sites  was  made  by  including  these 
new  absolute  measurements  in  Europe  to  the  IGSN  7 1 net.  The 
original  variance-covariance  matrix  of  IGSN  7 1 was  obtained  by 
single  precision  computations;  therefore,  no  more  than  six  new  sites 
were  obtained.  It  is  our  opinion  that  we  started  to  loose  signi- 
ficant figures  to  that  extent  that  the  seventh  selection  could  have 
been  effected  too  much  by  rounding  off  noise. 

We  repeated  the  above  mentioned  selection  processes  using 
the  UAU-net,  which  included  also  second  order  correction  terms 
to  the  calibrations  of  the  gravimeters.  The  second  order  effect  is 
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clearly  seen  from  the  preferred  selections  as  compared  to  the 
selections  obtained  from  IGSN  7 1 net. 


From  these  studies  it  became  clear  that  the  strong  correla- 
tions between  gravity  values  in  the  nets  were  disappearing  and  the 
existing  gravity  ties  with  their  current  accuracies  could  not  transfer 
the  effect  of  new  absolute  measurements  far  in  the  global  nets. 
Therefore,  in  the  future,  there  should  be  new  absolute  measurements 
of  gravity  on  different  continents  and  local  areas.  We  should  not 
try  to  make  new  global  adjustments  of  networks  but  rather  continental 
or  national  network  adjustments  including  new  absolute  measurements 
in  the  area  in  question. 

We  had  not  planned  originally  to  do  separate  analyses  for 
Africa  or  North  America  and  the  United  States.  However,  we  made 
extensive  study  in  the  areas,  especially  in  North  America  and  the 
United  States.  The  results  were  given  in  Chapter  6,  section  6.3. 

The  comments  were  given  in  the  same  section.  If  we  wish  to 
improve  the  U.  S.  Network,  we  must  do  more  accurate  measure- 
ments of  gravity  differences  between  the  stations  or  we  have  to 
make  a very  large  number  of  absolute  measurements  of  gravity. 

In  order  to  improve  accuracies  of  relative  gravity  measurements 
we  must  establish  a good  calibration  line  for  gravimeters  and 
possibly  improve  measuring  techniques  and  analyze  environmental 
effects  more  precisely  than  before. 

5.  Analyses  on  the  location  of  absolute  gravity  measurements  for 
calibration  of  gravimeters 

5.  1 Introduction 

During  the  adjustment  of  IGSN  71  net  it  was  recognized  that 
the  second  order  correction  term  to  the  calibrations  of  gravimeters 
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The  mathematical  model  for  inclusion  of 


r 


would  be  appropriate, 
this  second  order  term  was  worked  out  and  used  in  an  adjustment 
of  IGSN  7 1 net  (Uotila,  1974).  It  was  found  that  the  distribution 
of  absolute  sites  in  IGSN  7 1 was  not  good  for  solving  the  second 
order  term.  The  task  under  this  contract  was  to  determine  proper 
intervals  for  additional  absolute  sites  in  order  to  solve  the  second 
and  higher  order  terms  for  calibration  of  gravimeters. 

5.2  Mathematical  model 

A proper  mathematical  model  for  the  inclusion  of  second  or 
higher  order  correction  terms  to  the  calibration  of  gravimeters, 
must  include  the  dial  readings  because  the  higher  order  terms 
are  affected  by  the  location  of  the  readings  in  the  total  range  of  readings. 

A possible  mathematical  model  for  two  dial  readings  with  a 
gravimeter  in  a trip  including  the  third  order  correction  term  is: 

d»  - d»  + k»  (tt  - tj  ) + X®  (d}  - d*  ) + m*  (d*2  - dj2) 

+ n*  (d»  3 - d®  3)  - (g®  - g}  ) = 0 

where 

dj  , d(  = dial  readings  in  mgal  at  the  stations  i and  j, 

respectively,  corrected  for  all  known  systematic 
effects. 

k = coefficient  for  drift. 

t{ , t = time  of  observation  of  the  dial  readings  at  station 
i and  j,  respectively. 

X - coefficient  for  a linear  scale  factor  term, 
m = coefficient  for  a second  order  scale  factor  term, 
n = coefficient  for  a third  order  scale  factor  term. 
g(,  gj  = gravity  values  in  mgal  at  the  stations  i and  j, 
respectively. 

a = superscript  indicating  theoretical  or  adjusted  value. 
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We  would  have  a similar  equation  for  each  gravity  differ- 
ence. The  general  form  of  this  mathematical  model  is 

F (X®  , L»  ) = 0 

where  X*  = theoretical  or  adjusted  values  of  parameters 

La  = theoretical  or  adjusted  values  of  quantities  that  have 
been  observed  or  to  be  observed. 

In  this  case  Xa  would  include  k»  , iH,  m*  and  n»  for  each 
instrument  for  the  time  periods  during  which  they  are  considered 
invariant,  and  g?  for  each  station  included  in  the  net.  1 he  vector 
L*  would  have  each  dial  reading  of  the  instrument  as  an  element. 
The  values  t,  and  t are  considered  errorless  in  this  model.  The 
model  can  be,  of  course,  expanded  to  include  even  higher  order 
correction  terms  to  the  calibration  of  the  gravimeters. 


is: 


The  usual  minimum  variance  solution  for  the  above  model 
X = -[AT(BP'1BT)'*A]'1  A r (BP  *BT)  1 W 


where 


A = 

P'1 


— I ; B = — 1 
^X*  X*=X0'  *L,L'--Lb 


= £ - variance-covariance 


W F(Lb,X°) 

matrix  of  observed  quantities. 


Lb  = values  of  observed  quantities, 

X°  = app roxima te  values  of  parameters, 

X*  = X°  + X = adjusted  values  of  parameters. 

The  va  rianc e-cova rianc  e matrix  of  the  parameters  has  the  form: 


2^,  = [AT(BP'1BT)'lA]-‘: 

For  each  absolute  gravity  measurement  we  would  have  a mathemati- 
cal model: 

c*  - gt  = 0 
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where  c»  is  the  adjusted  value  of  the  absolute  gravity  measurement 
at  i station  and  g»  is  the  adjusted  gravity  value  of  the  i*”^1  station. 
The  absolute  values  can  be  added  to  the  earlier  expression: 

X = -[AT(BP  1 Bt) _1  A + Pj'1  [at(bp  ’b7)'1  w-pxw.] 

where  dimensions  of  Py  are  the  same  as  AT(BP'lBT)_,A  or  u x u 
when  the  dimensions  of  X are  u x 1.  All  other  elements  of  Py  are 
zero  except  those  diagonal  elements,  which  correpond  to  corrections 
to  the  g values  where  absolute  measurements  have  been  made. 

The  non-zero  diagonal  element,  px  , is  equal  to — -r—  which  is  the 

1!  CK  th 

reciprocal  of  the  variance  of  the  absolute  measurement  at  the  i 
station.  The  corresponding  wc  element  is 

w=t  = c*  - 8? 

where  c*  is  the  measured  absolute  value  at  the  ith  station  and  g°  is  the 

th 

approximate  gravity  value  at  the  i station  adopted  for  the  solution. 
All  other  elements  of  W? -matrix  are  zero  except  those  correspond- 
ing to  the  sites  where  absolute  measurements  have  been  made.  The 
corresponding  variance-cova  riance  matrix  of  parameters  is: 

rXt  = [at(bp~‘bt) "*  A + Pj  ‘ 

If  the  higher  order  drift  terms  or  some  environmental  factors,  such 
as  temperature  and  pressure  were  to  be  included  in  the  model,  they 
could  be  easily  added  by  modifying  the  mathematical  model  corre- 
spondingly. 

5.3  Results  of  analyses 
5.3.  1 Linear  correction  term 

If  we  were  interested  in  solving  only  the  linear  correction  term  to 
the  calibration  of  a gravimeter,  then  with  current  accuracies  of  absolute 


gravity  measurements,  10-20  ygal,  it  would  be  satisfactory  to  have 
only  two  absolute  measurements,  located  at  the  stations  having  the 
minimum  and  maximum  gravity  values  of  the  calibration  line.  The 
additional  information  obtained  by  including  more  absolute  measure- 
ments between  these  points  would  not  increase  the  accuracy  of 
determination  of  the  linear  calibration  correction  term  as  much  as 
making  these  additional  measurements  at  these  points  equally  divided 
between  the  two  points  of  the  calibration  line. 

5.3.2  Linear  and  the  second  order  correction  term 


In  the  case  that  we  wish  to  determine  the  linear  and  the 
second  order  correction  terms  to  the  calibration  of  the  gravimeters 
the  analyses  showed  that  we  should  have  absolute  measurements  at 
the  following  locations  along  the  calibration  line:  at  the  station 
having  maximum  gravity  value  and  at  the  station  having  minimum 
gravity  value  and  the  third  one  at  the  station  having  close  to  the 
average  of  minimum  and  maximum  gravity  values  along  the  calibration 
line.  Several  alternative  situations  were  examined,  but  this  simple 
system  seemed  to  give  the  best  results  for  the  case.  The  additional 
information  obtained  by  including  more  absolute  measurements 
between  these  points  would  not  increase  the  accuracy  of  the  deter- 
mination of  the  correction  terms  to  gravimeters  as  much  as  making 
these  new  measurements  at  these  three  points  or  at  the  vicinity  of 
these  three  points. 

5.  3.  3 Higher  order  correction  terms 

We  had  available  the  factory  calibration  tables  for  26 
La  Coste-Romberg  G gravimeters.  The  curves  of  the  original  cal- 
ibrations were  plotted  for  the  full  operational  scale  of  the  gravi- 
meters. In  table  12  we  see  the  results  of  the  fits  of  2nd,  3rd, 
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4th  and  5th  order  polynomials  to  the  calibration  carves.  Figures 
5-  12  give  some  sample  residuals  of  the  polynomial  fits  of  typical 
calibration  curves.  It  should  be  noted  that  the  scale  in  the  vertical 
axis  is  not  the  same  in  all  figures.  The  percentages  accounted 
by  the  polynomials  are  given  in  table  12.  The  percentage 
accounted  for  is  computed  using  the  following  formula 

Percentage  accounted  for  = Xl  ^ vi  x 100 


where  x1  = residuals  after  a first  order  polynomial  fit. 

Vj  = residuals  after  a particular  polynomial  fit. 

Similar  analyses  were  made  for  the  calibration  curves  for 
the  measuring  interval  which  has  been  used  in  IGSN  71  net  and  for 
the  interval  used  in  the  U.S.A.  The  corresponding  graphs  of  resi- 
duals were  also  plotted  for  these  polynomials,  but  in  much  larger 
scale.  Examples  are  given  for  two  gravimeters,  namely  L045  and 
L803  in  Figures  13-16.  We  should  again  note  differences  in  scale 
from  figure  to  figure.  We  tried  to  do  also  some  spectral  analyses, 
but  did  not  find  them  helpful  in  these  evaluations  at  this  time. 

Based  solely  on  the  examination  of  the  factory  calibration 
curves  we  can  conclude  that  the  curves  can  be  reproduced  to  the 
operational  accuracy  of  La  Coste-Romberg  G gravimeters  by  making 
absolute  measurements  at  500  mgal  intervals.  It  is  known  that 
these  gravimeters  have  been  calibrated  in  the  factory  using  200  mgal 
rider  at  various  parts  of  the  meter's  range  (Harrison  and  La  Coste, 
1978).  This  200  mgal  interval  does  not  seem  to  show  up  in  the 
calibration  curves.  It  is  clear,  of  course,  that  if  absolute  measure- 
ments are  made  at  200  mgal  interval  along  a calibration  line,  we 
should  be  able  to  control  calibration  of  the  gravimeters  to  the  same 
accuracy  as  the  factory  calibration  does;  however,  we  must  take 
into  consideration  changes  in  local  environmental  conditions,  such 
as  the  effects  of  tidal  variations,  changes  in  water  level,  etc. 


Table  12 

Analyses  of  LaCo ste- Rombe rg-fa cto ry  calibration  curves 


Gravimeter  Trend  Removed-Percent  Accounted  Order  of  Polynomial 


# 

2nd 

3rd 

4 th 

5 th 

L 00  1 

26.  6 

88.9 

99.  16 

99.  36 

L 002 

84.9 

9 1.  7 

97.  31 

99.  56 

L 007 

39.  7 

92.  0 

99.8  1 

99.85 

L 009 

52.  4 

97.  5 

99.95 

99.  95 

L 011 

15.  7 

83.  9 

99.  74 

99.  75 

L 0 12 

2.  5 

94.9 

99.93 

99-94 

L 020 

83.8 

96.  2 

99.97 

99.98 

L 043 

94.  5 

97.8 

99.9  1 

99.9  1 

L 044 

84.  9 

92.  5 

99.72 

99.80 

L 045 

19.  7 

90.  7 

95.  08 

96.  74 

L 046 

47.  1 

94.  1 

99.  37 

99.43 

L 047 

16.4 

60.  4 

99.  58 

99.82 

L 048 

85.  6 

89.  3 

99.  78 

99.97 

L 050 

9 1.9 

92.  0 

99.71 

99.94 

L 056 

28.2 

64.  1 

99.00 

99.69 

L 057 

79.  7 

94.  7 

99.  70 

99.90 

L 074 

65.2 

90.  2 

99.94 

99.95 

L 075 

87. 2 

90.  2 

99.87 

99.87 

L 093 

86.4 

98.6 

99.99 

99.99 

L 115 

63.  0 

88.  7 

99-94 

99.99 

L 122 

85.4 

96.0 

99-98 

99.  98 

L 137 

0.6 

69.  3 

99-  29 

99.70 

L 140 

9 1.8 

96.9 

99.89 

99.98 

L 803 

97.  8 

98.  2 

99.  34 

99.45 

L 808 

62.  0 

9 1.  3 

99.  77 

99.99 

L 903 

93.  7 

95.  0 

99.  58 

99.62 
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(Brein,  et  al,  1977).  The  periodic  errors  in  the  measuring  screw 
should  be  also  modeled  if  at  all  possible.  These  errors  are  caused 
by  excentricity  in  the  screw  resulting  in  wobble  and  non-linearity  of 
lever  system.  The  screw  problems  result  in  errors  typically  35 
^gal  in  amplitude  in  the  G-meters  once  or  twice  per  turn  of  the 
screw  (about  70  and  35  mgal)  (Harrison  and  La  Coste,  1978). 

What  type  of  calibration  lines  and  how  many  absolute  measurements 
are  required  to  control  this  kind  of  behavior  are  the  questions,  for 
which  the  answers  are  outside  of  the  scope  of  this  research,  but 
should  be  studied  in  the  future.  It  has  been  suggested  (Brein,  et 
al,  1977)  that  there  should  be  1 mgal  baselines  with  5 to  10  inter- 
mediate stations  for  this  purpose.  How  many  of  this  kind  of  baselines 
are  needed  and  do  we  need  absolute  measurements  in  these,  should 
be  also  examined  in  the  future.  When  other  types  of  gravimeters 
such  as  La  Coste-Romberg  D-meter  are  considered,  we  must 
recognize  that  the  gravity  ranges  differ  from  those  of  G-meter. 

The  D-meter  measures  only  about  a 200  mgal  change. 

5.4  Conclusions 


If  we  were  to  control  only  a linear  correction  term  to  the 
calibration  of  a gravimeter  and  accept  the  factory  calibration  for 
the  local  variations,  two  locations  of  the  absolute  measurements 
would  be  the  most  favorable  ones,  namely  the  stations  having 
minimum  and  maximum  gravity  values  among  the  stations  in  the 
calibration  line.  If  we  were  to  control  also  the  correction  to  the 
second  order  term,  an  additional  station  should  be  selected,  which 
would  have  a gravity  value  half  way  between  maximum  and  minimum 
gravity  values  along  the  calibration  line.  If  we  wish,  for  some 
other  reasons,  to  select  more  than  these  stations,  they  should  be 
equally  divided  in  the  vicinities  of  the  selected  ones.  The  "vicinity" 
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means  in  this  context  - having  close  to  the  same  gravity  values. 


If  we  wish  to  reproduce  similar  calibration  curve  as  the  one 
supplied  by  the  factory  for  G-meters,  S00  mgal  interval  is  satis- 
factory provided  that  environmental  effects  are  taken  into  consider- 
ation. If  we  wish  to  obtain  10  pgal  or  better  accuracy  in  measure- 
ments of  gravity  differences,  we  must  establish  calibration  lines  of 
smaller  gravity  differences  related  to  the  gravity  differences  to  be 
measured.  For  accurate  measurements  and  field  calibrations 
continuous  recording  of  tidal  and  other  variations  of  gravity  should 
be  made  at  the  stations  for  correcting  values  to  the  normal  values. 
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